DUDLEY KITOX LIERARY 
UAVAL PO I'L. r , •“ Y'T ■■’CHOOL 
ItOilTLBEY, o.uTxii. 9SG4S^&^ 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 



DEVELOPMENT OF A MATHEMATICAL MODEL 
THAT SIMULATES THE LONGITUDINAL, AND 
LATERAL-DIRECTIONAL RESPONSE OF THE 
F/A-18 FOR THE STUDY OF 
FLIGHT CONTROL RECONFIGURATION 

by 

Fredric W. Rojek 
September 1 986 



Thesis Advisor: Daniel J. Collins 



Approved for public release; distribution unlimited. 



J232476 



:URirY CLASSiFiCATIQN OF THIS PAGE 



REPORT DOCUMENTATION PAGE 



REPORT ScCLRlTY CLAS5iFiCAT!ON 

JNCLASSIFIED 



lb. RESTRICTIVE MARKINGS 



SECURITY CLASSiF'CATION AUTHORITY 



DECLASSIFICATION / DOWNGRADING SCHEDULE 



3 DISTRIBUTION /AVAILABILITY OF REPORT 

Approved for public release; distri- 
bution unlimited. 



PERFORiM.’NG ORGANIZATION REPORT NUM3ER{S) 



5. MONITORING ORGANIZATION REPORT NUM3ER{S) 



NAME OF PERFORMING ORGANIZATION | 

^aval Postgraduate School! 



6b. OFFICE SYMBOL 
(If applicable) 

Code 67 



7a. NAME OP MONITORING ORGANIZATION 

Naval Postgraduate School 



ADDRESS (City, State, and ZIP Cade) 
lonterey, California 



93943-5000 



7b. ADDRESS (Oty, State, and ZIP Code) 
Monterey, California 



93943-5000 



•NAME OF FUNDING /SPONSORING 
ORGANIZATION 



8b OFFICE SYMBOL 
(If applicable) 



9. PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 



ADDRESS (City, State, and ZIP Code) 



10 SOURCE OF FUNDING NUMBERS 



PROGRAM 


PROJECT 


TASK 


WORK UNIT i 


ELEMENT NO 


NO 


NO 


ACCESSION NO | 



TITLE f/nc/ude Secunty C/asi/fef/on; DEVELOPMENT OF A MATHEMATICAL MODEL THAT SIMULATES THE 
[)NGITUDINAL, AND LATERAL-DIRECTIONAL RESPONSE OF THE F/A-18 FOR THE STUDY 
FLIGHT CONTROL RECONFIGURATION 



personal AUTHOR(S) 
ROJEK, FREDRIC W. 



ji TYPE OF REPORT 


13b TIME COVERED 


14 DATE OF REPORT (Year, Month. Day) 


15 PAGE COUNT 




lister’s Thesis 


FROM TO 


September 1 986 


286 


i 

3 



SUPPLEMENTARY NOTATION 



COSATI COOES 



FIELD 



GROUP 



SUB-GROUP 



18. SUBJECT TERMS (Continue 0.1 reverse if necessary and identify by block number) 

F/A-18, Simulation, Reconf igurable , Flight, 
Controls 



\ ABSTRACT (Continue on reverse if necessary and identify by block number) j 

A linearized mathematical model is developed which simulates the dynamic | 
iBsponse of the Navy F/A-18 for the study of flight control reconfiguration. ] 
jie aircraft is modeled as a multi-input multi-output, sampled data, closed 1 
ystem, which couples the dynamics of the flight control system to the air- \ 
jraft linearized small perturbation equations. The discrete time, state j 
iriable equations for the system are then formulated. A computer program is.j 
pveloped which will compose the model matricies and compute the response of 3 
jie aircraft to stick and rudder inputs. | 
I To study flight control reconfiguration, the model allows individual I 
ptuation of either a left or right control surface. Aircraft response to | 
jie actuation loss of either the left or right stabilator is simulated in the| 
jrogram. The program is designed to implement the reconf igurable control | 
*^sleA under study for the Self-Repairing Digital Flight Control | 



1 distribution /AVAILABILITY OF abstract 

SXUNCLASSIFIED/UNLIMITED □ SAME AS RFT □ DTIC USERS 


21. ABSTRACT SECURITY CLASSIFICATION | 

UNCLASSIFIED | 


S NAME OF RESPONSIBLE INDIVIDUAL 

aniel J. Collins 


22b TELEPHONE (/nc/uG(e Area Code) 

(408) 646-2826 


22c. OFFICE SYMBOL 1 

Code 6 7 Co | 



I) FORM 1473, 84 MAR 



33 APR edition may be used until exhausted 
All other editions are obsolete 
1 



SECURITY CLASSIFICATION OF THiS PAGE 



SECURITY CLASSIFICATION OF THIS PAGE Daim Enf«f»d) 



Block 19. ABSTRACT (cont'd) 

The computer simulation was written in VS FORTRAN. A copy 
of the program and simulation results are included in the 
appendices . 



s N 0)02- LF- 014- 6601 

2 



security classification of this PAGEfW7i»n Data Bnfrtd) 



Approved for public release; distribution unlimited. 

Development of a Mathematical Model that Simulates the 
Longitudinal, and La teral-Di rect iona J Response of the F/A-IP 
for the Study of Flight Control Reconfiguration 



by 



Fredric W. Rojek 
Lieutenant, Unitec^'^States Navy 
B.S.E.E., State University of New York at Buffalo 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 



from the 

NAVAL POSTGRADUATE SCHOOL 
September 1 986 



ABSTRACT 



A linearized mathematical model is developed which simu- 
lates the dynamic response of the Navy F/A-18 for the study 
of flight control reconfiguration. The aircraft is modeled 
as a multi-input mu 1 t i -o u t p u t , sampled data, closed system, 
which couples the dynamics of the flight control system to 
the aircraft linearized small perturbation equations. The 
discrete time, state variable equations for the system are 
then formulated. A computer program is developed which will 
compose the model matricies and compute the response of the 
aircraft to stick and rudder inputs. 

To study flight control reconfiguration, the model allows 
individual actuation of either a left or right control sur- 
face. Aircraft response to the actuation loss of either the 
left or right stabilator is simulated in the program. The 
program is designed to implement the r econf igur able control 
mixer, currently under study for the S e 1 f - Re pa i r i ng Digital 
Flight Control System. 

The computer simulation was written in VS FORTRAN. A 
copy of the program and simulation results are included in 
the appendices . 
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INTRODUCTION 



I . 

Tactical aircraft face airborne and ground based threats 
which continue to grow in number and capability. In future 
combat engagements the lethality of these systems will pose 
a considerable threat to aircraft survivability. In addition 
to combat losses, the loss of aircraft in battle damage re- 
pair, or awaiting repair, will significantly reduce our tac- 
tical forces. Projections on the survivability of NATO 
forces [Ref. 1] during the initial days of engagement indi- 
cate 68% of the tactical aircraft will be out of action 
after the third day of battle. Of this total, 22% will be 
lost in combat and 46% lost while in repair or awaiting re- 
pair. This is indicated in Fig. 1.1 which comes from Refer- 
ence 1. It is clear from projections such as this, if our 
forces are to remain a superior threat to the enemy, contin- 
ued emphasis must be given to reduce the combat vulnerability, 
and increase the reliability and maintainability of our tac- 
tical aircraft. 

To improve the combat effectiveness of future tactical 
aircraft, the United States Air Force initiated the Self 
Repairing Flight Control System, Reliability and Maintain- 
ability Program. Reference 2 outlines the program plan and 
goals. An Air Force sponsored study [Ref. 1] showed that 
significant improvements in aircraft survivability and 
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reliability can be achieved with a self repairing digital 
flight control system. 

One of the principle features of a self repairing system 
is flight control reconfiguration. A reconf igurable flight 
control system automatically counters aircraft loss of con- 
trol due to impairment, or loss, of a control effector. The 
fundamental idea is to utilize the existing undamaged effec- 
tors to preserve the normal flying qualities of the unimpaired 
aircraft. The approach outlined in Reference 1 is to design 
a reconf igurable control mixer to be placed between the flight 
control laws and the control effectors. In the face of air- 
craft impairment the mixer would reallocate control commands 
to the unimpaired effectors so that flight critical pitch, 
yaw, and rolling moments would be preserved to the maximum 
extent possible. Using the control mixer concept, the exist- 
ing flight control laws would not be altered. 

Although the development of the self repairing system is 
intended for the advanced tactical aircraft, it is possible 
that the mixer could be implemented in existing airframes. 

This could be done with a control reconfiguration module 
interposed between the aircraft flight control computers and 
fly-by-wire actuators [Ref. 1]. 

For this thesis the McDonnell F/A-18 was chosen to study 
the reconf igurable control law concept utilizing the control 
mixer as described above. This choice was made based on the 
following assumptions given in Reference 1 for implementation 



\ 
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of the reconfiguration control law process in an existing 
airframe : 

1) ”The effector complement provides redundant effector 
systems and surplus control power for each flight 
critical control force and moment.” 

2) ”The flight control is a full authority, fly-by-wire 
digital flight control system.” 

3) ”A control law design exists which has been carried 
out for the unimpaired airplane, but which is suf- 
ficiently robust that only first order impairment- 
induced stability derivative changes need be 
accounted for in a drop-in r e co n f i g u r a b 1 e mixer.” 

Based upon these assumptions it is felt that the F/A-18 is 

well suited for the reconfigurable control law study. 

In this thesis a linearized mathematical model was de- 
veloped which simulates the dynamic response of the F/A-18 
aircraft to stick and rudder inputs. The model includes the 
control laws for the longitudinal, lateral, and directional 
axes for the cruise phase of flight. The fly-by-wire actua- 
tors and sensor dynamics were modeled and joined to the air- 
frame linearized, small perturbation model. The perturbation 
model was obtained from the Simulation Control and Technology 
Group, Flight Systems Branch, Strike Aircraft Test Direc- 
torate, at the Naval Air Test Center. A computer program 
was developed which, with existing software at the Naval 
Postgraduate School, will simulate F/A-18 longitudinal, lat- 
eral, and directional response to stick and rudder inputs. 

The model was designed to implement the control mixer 
gain concept outlined above. Future work at NPS will utilize 
the simulation program to develop algorithms for determining 
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the control mixer gain matrix. In addition the program will 
be used to study modern control augmentation systenrs and air- 
craft stability and control. 
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II. MODEL METHODOLOGY 



A. INTRODUCTION 

This chapter discusses the methods and assumptions used 
to formulate a mathematical model which simulates the dynam- 
ic response of the F/A-18 aircraft. The flight control sys- 
tem is described, including the simplifying assumptions used 
to develop the control system model. An overview of the com- 
plete system, which couples the control system model with 
the airframe small perturbation model, is then given with a 
brief description of each functional component. 

B. FLIGHT CONTROL SYSTEM DESCRIPTION 

A detailed description of the F/A-18 control system and 
theory of operation can be found in the Flight Control Sys- 
tem Design Report by McDonnell Aircraft Company [Ref. 3]. 

The following discussion briefly describes the basics of the 
flight control system and the control law mechanization, and 
is intended to facilitate understanding the model development. 

The primary flight control system in the F/A-18 is a 
fly-by-wire, full authority, control augmentation system. 

The control law computations are performed by four flight 
control computers operating in parallel. Each computer re- 
ceives input from the aircraft motion sensors, air data 
computer, and pilot stick commands. The computer operates 
on the input signals according to the control law algorithms 
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and outputs the command signals to fly-by-wire electrohy- 
draulic ser voac tua tor s . Figure 2.1, taken from Reference 3, 
shows a functional block diagram of the flight control system. 
Exclusive of angle of attack and air data sensors, the system 
has quadraplex redundancy. The system provides two fail op- 
erate performance for augmented motion feedback control. A 
third failure causes the system to revert to either open loop 
direct electrical link control, or stabilator mechanical 
control . 

The control augmentation system is gain scheduled with 
angle of attack and air data to provide optimum flying quali- 
ties throughout the flight envelope. Cross axis intercon- 
nects (e.g., rolling surface to rudder interconnect) are 
provided for turn coordination and maneuverability at high 
angles of attack. The control system also provides feedback 
to counter inertial coupling at high roll rates. 

The F/A-18 has ten primary flight control surfaces: 

Right and left stabilators, leading edge flaps, trailing 
edge flaps, ailerons, and rudders. Longitudinal control is 
provided by collective stabilator, and collective leading 
and trailing edge flaps. La t e r a 1 -d i r ec t io na 1 control is 
provided by differential stabilator, differential leading 
and trailing edge flaps, ailerons, and rudders. Collective 
leading and trailing edge flap deflections are scheduled by 
the control laws and are a function of angle of attack. The 
flap positions are designed to provide optimum L/D during 
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Ps = Static pressure 

Qc = Dynamic pressure 



© 



= Trim motor 



— » » » » = Indicates level of electrical redundancy 



Figure 2.1 Functional Block Diagram of Flight Control System 
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cruise, and improve flying characteristics in maneuvering 
and high angle of attack flight. Figure 2.2, copied from 
Reference 3, shows the control surface positions and direc- 
tion for positive deflection. 

The objective of the thesis did not require a complete 
model of the F/A-18 control system as given in Reference 3. 
The following assumptions were made to reduce the model 
complexity . 

1) The aircraft is operating in the up and away flight 
phase. Under this assumption the control laws are 
operating in the a u t o- f 1 a p s-u p configuration. Con- 
trol law configurations for the takeoff and landing 
phases were not modeled. 

2) Only inner loop control is modeled. For inner loop 
control the pilot provides commands to the system. 

The auto functions (outer loop control) were not 
modeled . 

3) Control is provided by the control augmentation system. 
The unaugmented modes such as direct electrical link 
or mechanical backup are not modeled. 

4) The failure logic provided to reconfigure the control 
laws in the event of a sensor or actuator failure is 
not mod e 1 ed . 

5) The aircraft is operating with gear up, speedbrakes 
in, and no external stores. 

6) Spin mode control logic is not modeled. 

7) The aircraft trim system is not modeled. 

8) High angle of attack conditions are not considered in 
the thesis model, therefore control law configurations 
for this flight condition are not modeled. For this 
thesis high angle of attack is defined as flight con- 
ditions above 15 degrees. 
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C. F/A-18 DYNAMIC MODEL OVERVIEW 

The functional block diagram of the F/A-18 model which 
couples the flight control system to the basic airframe is 
shown in Figure 2.3. The d ia g r am r e p r e s en t s a multi-input 
multi-output, sampled data, closed loop control system. 

Theory on the analysis of sampled data systems, as the one 
shown in Fig. 2.3, is extensive and covered in a number of 
texts (see references). The development of the F/A-18 model 
assumes the reader has a rudimentary understanding of control 
theory, and in particular the theory of sampled data control 
systems . 

In the nomenclature used to represent the control signals 
pilot inputs are prefixed by the letter ’P'. Actuating sig- 
nals, and signals produced by the aircraft sensors are pre- 
fixed by the letter ’E'. Control surface deflections are 
prefixed by the letter 'D'. The nomenclature used to repre- 
sent the aircraft perturbed motion variables and control 
surfaces is given in Table 2.1. To denote the motion axis 
which is being controlled, the signal \\rill be suffixed with 
either x, y, or z to denote longitudinal, lateral, or direc- 
tional axis respectfully. Finally a matrix or a vector will 
be denoted by an upper case letter. A scalar will be denoted 
by a lower case letter. 

The input vector shown in Fig. 2.3, 

= [px(t) py(t) pz(t)] (2.1) 
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TABLE 2.1 



MOTION VARIABLES AND CONTROL SURFACE NOMENCLATURE 



q 


pitch rate 


DEG/SEC 


nz 


normal acceleration 


G 


aa 


angle of at tack 


DEG 




yaw rate 


DEG/SEC 


r r 


roll rate 


DEG/SEC 


ny 


lateral acceleration 


G 


s tr 


right stabi la tor 


DEG 


s 1 1 


left s tabilator 


rt 


ter 


right trailing edge flap 


?T 


tel 


left trailing edge flap 


If 


ler 


right leading edge flap 


IT 


lei 


left leading edge flap 


TT 


a r 


right aileron 


Tf 


al 


left aileron 


If 


r r 


right rudder 


IT 


rl 


left rudder 


If 



represents the aircraft longitudinal and lateral stick, and 
rudder deflection in inches. The output vector, 

Y(t) '" = [q(t) nz(t) aa(t) yr(t) rr(t) ny(t)] (2.2) 

represents the perturbed motion of the aircraft about some 
steady state operating condition. The motion variable units 
are degrees, degrees/sec, and G^s. Each block in the dia- 
gram contains a mathematical model which simulates the dynam- 
ics of that particular component. The control law block 
contains the aircraft flight control law algorithms modeled 
as linear, time invariant, discrete state equations. One 
processing channel of the flight control computer described 
above is represented in the control law model. The input 
vectors to the control law model are the discrete stick and 
rudder input signals in inches. 
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[px(k) py(k) pz(k)] 



(2.3) 



P(k) 



and the discrete motion feedback signals from the sensors, 

Es(k)*' = [q(k) nz(k) aa(k) yr(k) rr(k) ny(k)] (2.4) 



(The units of the mo t i o n f ee d back s i gna 1 s are degrees, degrees/ 
sec, and G^s). The input vector to the control law equations 
is therefore 

U(k)^ = [Es(k) I P(k)] (2.5) 

The output vector from the control law block, 

Ec(k)^ = [estr(k) estl(k) eler(k) elel(k) eter(k) 

etel(k) ear(k) eal(k) err(k) erl(k)] (2.6) 

represents the discrete command signals to the flight con- 
trol actuators in degrees. The actuator command signals en- 
ter the GAIN block which represents the configuration gain 
matrix discussed in the introduction. The state space equa- 
tions shown in the actuator block model the dynamics of the 
flight control actuators. The input vector to the actuator 
block , 

Ea(t)^ = [estr(t) estl(t) eler(t) elel(t) eter(t) 

etel(t) ear(t) eal(t) err(t) erl(t)] (2.7) 

represents the continuous time, actuator command signals in 
degrees. The output vector from the actuator block 
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represents the control surface deflections in degrees, 

D(t)^ = [dstr(t) dstl(t) dler(t) dlel(t) dter(t) 

dtel(t) dar(t) dal(t) drr(t) drl(t)] (2.8) 

The deflection vector is input to the airframe small pertur- 
bation model represented by the state variable equations in 
the airframe block. The output vector from the small per- 
turbation model, Y(t), enters the sensor dynamics block which 
contains the state variable model for the aircraft rate gyros, 
accelerometers, and angle of attack sensors. As discussed 
above the sensors output the feedback signals which are sent 
to the control laws. In Chapter III the mathematical models 
which simulate the dynamics of each component in Fig. 2.3 
will be developed in detail. 

Analog to digital converters are modeled as impulse sam- 
plers. It is assumed that all samplers are operating at the 
same, constant sampling rate. (The actual system uses multi- 
rate sampling. In the thesis model only a single sampling 
rate is used. The program actually allows any desired sam-- 
pling rate to be input.) The mathematical operation per- 
formed by the impulse samplers is shown in Figure 2.4. 

Modeling the analog to digital converters as impulse samplers 
is a valid assumption if the quantization error of the actual 
system is at an acceptable level. Digital to analog conver- 
ters are modeled as zero order hold devices. The mathemati- 
cal operation performed by the zero order hold is shown in 
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Fig. 2.5. Finally it is assumed that the processing time be- 
tween the sampler inputs and zero order hold output is very 
much less than the sampling period (i.e., no processing delay 
time is assumed). This assumption is used in transforming 
the continuous time state equations (actuators, airframe, and 
sensors) into discrete time equations. 
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Figure 2.4 Mathematical Operation Performed by Impulse Sampler 
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Mathematical Operation 



Performed by 



Zero Order Hold 



26 



III. MODEL DEVELOPMENT 



A. CONTROL LAW MODEL DEVELOPMENT 

The longitudinal, lateral, and directional control law 
models were developed from the block diagrams and information 
provided in Reference 3. In addition to the simplifying as- 
sumptions listed in Chapter II, many of the components in 
this system were eliminated based on the following 
considerations : 

1) The control law model will be coupled to the small 
perturbation model of the F/A-18. It is assumed 
that operation of the control system will remain 
within a linear region. Therefore the non-linear 
components in the control laws were eliminated. 

These include position limiters, rate limiters, 
and dead band regions. Other non-linear functions 
in the system which were essential to the model 
(pitch stick gradient for example) were linearized 
by a taylor series expansion. 

2) For the same reasons discussed above, the portions 
of the system which provide inertial decoupling 
were not modeled. 

3) In the model, the control signals are input as dis- 
crete signals. Therefore stick and rudder dynamics 
are not modeled. 

4) Noise, which may be detrimental to control system 
performance due to aliasing, is not introduced into 
the model. Therefore the anti-aliasing prefilters 
were not included in the model. 

5) The structural modes were not included in the F/A-18 
airframe model. Therefore the structural notch fil- 
ters were not included in the model. 

6) To prevent discontinuities in the signals, the con- 
trol laws utilize faders in portions of the system. 
Discontinuities could occur during start up, fail- 
ures, or transitions. Since none of these condi- 
tions are included in the thesis model, faders have 
not been modeled. 
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7) The F/A-18 control system uses multi-rate sampling 
in the input and feedback paths (20, 40 and 80 hz 
sampling rates are used). To develop a state var- 
iable model only a single sampling rate was con- 
sidered. Therefore the iteration averagers, used 
to mathematically combine two discrete signals of 
different sampling periods were not modeled. The 
simulation program allows any desired sampling rate 
to be input. For this thesis 80 hz was used as the 
sampl ing r a t e . 

Figures 3.1 and 3.2 show the simplified block diagrams 
of the longitudinal and lateral-directional control laws. 

The inputs are the discrete stick and rudder commands, and 
motion feedback signals. The output signals are the commands 
to be sent to the flight control actuators via the reconfig- 
urable gain matrix. Together Figs. 3.1 and 3.2 make up the 
control law block shown in Fig. 2.3. 

The blocks in the control law model represent two basic 
transfer functions: Function gains, and digital filters. 

Table 3.1 lists the notation used to represent the transfer 
functions. Function gains and digital filters are described 
below. 



. TABLE 3 . 1 



CONTROL LAW TRANSFER FUNCTION NOTATION 
Prefix Transfer function 



F Function gain 

P Longitudinal digital filter 

R Lateral digital filter 

Y Directional digital filter 
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Figure 3.1 Simpified Longitudinal Control Law Diagram 
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Functions 



1 . 

The functions in the control law diagrams perform the 
system gain scheduling described in Sec. II. B. The functions 
operate on the air data, angle of attack, and normal acceler- 
ation to compute the system gains. Table 3.2 lists the nota- 
tion used to represent the function inputs. The mathematical 
equations which define the gain schedules are given in Refer- 
ence 3. Copies of the functions used in the control law 
model are given in Appendix A. In the simulation program the 
function gains are computed using steady state conditions for 
all input values. 



ps 
qc 
r i 
n 2 

alpha 



TABLE 3 . 2 

FUNCTION INPUT VARIABLES 



Indicated static pressure lbs/ft2 
Dynamic pressure lbs/ft2 

Pressure ratio (ps /qc) ND 

Normal acceleration G 

Angle of attack DEG 



2 . Digital Filters 

Lead-lag filters are used in the system to shape the 
output response and provide adequate gain and phase margins. 

An integrator is used in the forward loop of the longitudinal- 
system to provide zero steady state error between command and 
feedback. The control system design report [Ref. 3] gives 
the filter^s continuous time transfer function. For the con- 
trol law model, the digital filter coefficients were computed 
using the Tustin transform. All filters were modeled as first 
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order systems with two numerator coefficients and a single 
denominator coefficient. For example 

H(Z) = P9N1^--Z + P9N2 (3.1) 

Z - P9D 

is the model for filter P9. Appendix B gives the digital 
filters used in- the model, and the method used to -compute the 
coefficients. 

3 . State Space Models 

Standard control system analysis techniques were used 
to derive the state space models for longitudinal and lateral- 
directional control laws shown in Figs. 3.1 and 3.2. The 
following sections outline the procedures used, 
a. Longitudinal System 

The block diagram in Fig. 3.1 contains five 
input/output signal flow paths: 

1) Pitch rate to collective stabilator 

2) No r ma 1 a c c e 1 e r a t i o n to collective stabilator 

3) Longitudinal stick to collective stabilator 

4) Angle of attack to collective leading edge flap 

5) Angle of attack to collective trailing edge flap 

To obtain the individual path transfer functions the signals 
are mathematically combined to give the following three Z- 
transform equations: 

Estx(Z) = Hl(Z)Eq(Z) - H2(Z)Enz(Z) - H3(Z)Epx(Z) (3.2) 
Elex(Z) = H4(Z)Eaa(Z) (3.3) 
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Etex(Z) = H5(Z)Eaa(Z) 



(3.4) 



The expressions for each transfer function are given in Ap- 
pendix C, To obtain a state space expression for longitudi- 
nal control laws, the individual Z-transfer functions are 
first expressed in state variable form. The state variable 
equations are then combined according to equations 3.2 - 3.4 
This procedure is as follows: 

1) The transfer function for the pitch rate to collective 
stabilator path can be expressed as: 

Hl(z) = bOZ^ + blZ + b2 (3.5) 

(Z-P9D) (Z-P2D) 



The numerator coefficients are functions of the system gains 
and filter coefficients. The roots in the denominator are 
the poles from filters P9 and P2. (Appendix C gives the de- 
tailed expressions for the numerator coefficients). 

2) The state space representation of Eq . 3.5 is obtained 

using the parallel programming method outlined in 
Reference 4 . 



‘xl(k + l)" 




■P9D 0 ■ 


"xl(k)' 


1 


'1“ 


_x 2 ( k + 2 )_ 




0 P2D 


_x2(k)_ 


+ 


_ 1 _ 



eq(k) 



(3.6) 



estxl(k) =|qstl qst 



2 ] 



xl(k)- 
_x 2 ( k ). 



q s 1 3 



eq(k) 



(3.7) 



Applying the same procedures to the remaining Z-transfer 
functions in Eq. 3.2 results in similar expressions: 
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"x3 ( k + 1 )” 




■P9D 0 ■ 


x3(k)- 




’I" 


_x4 ( k + 2 )_ 


= 


. 0 P5D 


_x4 ( k )_ 


+ 


_ 1 _ 



e n z ( k ) 



(3.8) 



estx2(k) = nzstl nzst2 



'xl(k) 

.x2(k)J 



+ nzst3 enz(k) (3.9) 



for the normal acceleration path, H2(Z); and 



x-5(k + l) = P9D x5(k) + 1 px(k) 
estx3(k) = pxstl x5(k) + pxst2 px(k) 



(3.10) 

(3.11) 



for the longitudinal stick path, H3(Z). 

Appendix C details the procedures used to com- 
pute the coefficients in the output equations. The nomencla- 
ture used to represent the coefficients in the output equation 
combines the notation of the input signal and output 
control surface, followed by a number indicating the coeffi- 
cient's numerical order in the equation. For example: 



qs tl , qs t2 , qs t3 



are the first, second, and third coefficients in the pitch 
rate to stabilator output equation (Eq. 3.7). With this 
system of nomenclature the respective signal path of the co- 
efficient is easily identified. 

3) The state variable equations for H1(Z), H2(Z) and H3(Z) 
are now combined according to Eq . 3.2. 

estx(k) = estxl(k) - estx2(k) - estx3(k) (3.12) 
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Adding Eqs. 3.7, 3.9 and 3.11 gives 



e s t X ( k ) 



= [qstl 



qst2-nzstl— nzst2 



- p X s 1 1] 



fx 1 ( k )- 
x2(k) 
X 3 ( k ) 
x4(k) 
J< 5 ( k 



+ 



[qst3-nzst3] 



eq(k) 
enz ( k ) 



+ “pxstl px(k) 



(3.13) 



for the output equation. Equations 3,6, 3.8, and 3.10 can 

be combined to give 
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eq(k) 
enz ( k ) 
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o 
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. 0 . 



px (k) 

(3.14) 



The state space representation of equation 3.2 is given by 
Eqs. 3.13 and 3.14. These equations output the collective 
stabilator command given the longitudinal stick and motion 
feedback inputs. Note that the motion feedback and longitu- 
dinal inputs have been separated. This facilitates coupling 
the control law equations to the aircraft equations to be 
developed later . 

4) Similar state variable equations are derived for the 
angle of attack to collective flap path transfer 
functions H4(Z), and H5(Z): 

x6(k + l) = PHD x6(k) + 1 eaa(k) (3,15) 

elex(k) = aalel x6(k) + aale2 eaa(k) (3.16) 



for the AOA to collective leading edge path H4(Z); and 
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x7(k+l) = P12D x7(k) + 1 eaa(k) 
etex(k) = aatel x7(k) + aate2 eaa(k) 



(3.17) 

(3.18) 



for the AOA to collective trailing edge path H5(Z). 

5) The state equations for collective stabilator, and 
collective leading and trailing edge flaps are now 
combined to give the following state variable model 
for the longitudinal control laws: 
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(3.19) 
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(3.20) 



Equations 3.19 and 3.20 are the state space representation 
of the longitudinal control laws shown in the block diagram 
in Fig. 3.1. The equations represent the longitudinal con- 
trol law model which computes the collective stabilator 
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command, and collective flap commands, given the longitudinal 
stick and mo t ion feedback inputs. 

b. Lateral-Directional System 

The procedures outlined above are applied to the 
Fig. 3.2 to obtain the state space model for the lateral- 
directional system. (To complete the discussion of the con- 

V 

trol law model these procedures will be briefly described.) 

The individual path transfer functions are first obtained by 
mathematically combining the signals in Fig. 3.2 to give the 
following equations: 

Esty(Z) = -H6(Z)Err(Z) + H7(Z)Py(Z) + H8 ( Z ) Pz ( Z ) ( 3 . 21 ) 
Eley(Z) = -H9(Z)Err(Z) + H10(Z)Py(z) (3.22) 

Etey(Z) = -HlKZ)Err(Z) + H12(Z)Py(z) (3.23) 

Ea(z) = -H13(Z)Err(Z) + H14(Z)Py(Z) + H15(Z)Pz(Z) (3.24) 

Er(Z) =-Hl6(Z)Eyr(Z) t H17(Z)Err(Z) + 

H18(Z)Eny(Z) + H19(Z)Py(Z) + H20(Z)Pz(Z) (3.25) 

These transfer functions represent the following input/output 
signal, paths, numbered respectfully: 

6) Roll rate to differential stabilator 

7) Lateral stick to differential stabilator 

8) Rudder pedal to differential stabilator 

9) Roll rate to differential leading edge flap 

10) Lateral stick to differential leading edge flap 

11) Roll rate to differential trailing edge flap 

12) Lateral stick to differential trailing edge flap 
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13) 



Roll rate to aileron 



14) Lateral stick to aileron 

15) Rudder pedal to aileron 

16) Yaw rate to rudder 

17) Roll rate to rudder 

18) Lateral acceleration to rudder 

19) Lateral stick to rudder 

20) Rudder pedal to rudder 



The expression for each transfer function and corresponding 
state equation are given in Appendix C. Note that the dif- 
ferential stabilator, ailerons, and rudder contain the trans- 
fer functions for the cross axis interconnects (e.g., H8(Z) 

and H15(Z) represent the rudder to rolling surface intercon- 
nect, and H19(Z) represents the rolling surface to rudder 
interconnect, ) 

The state equations for the lateral and direction- 
al control laws are given as: 
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(3.26) 



for the lateral system; and 
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(3.27) 
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for the directional system. Note that the lateral system is 
of order zero. No filters were included in the model, 
c. Combining the Control Law Models 

Equations 3.19 and 3.20, and 3.26 - 3.28 are com- 
bined to give the 3-axis control law model. 
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Equations 3.29 and 3.30 are written in terms of the matrix 
coefficients in Appendix D. These equations are represented 
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by the discrete state equations in the control law block in 
Fig. 2.3. 

3 . Variable Gain Matrix 

As a final step in the control law model dev^elopment, 
the variable gain matrix is introduced. The command signals 
which are output from the control laws (Eqns. 3.29 and 3.30) 
are distributed to the right and left actuators according to 
the diagram in Fig. 3.3. For the unimpaired aircraft the 
individual gains in Fig. 3.3 are set to unity. • The follow- 
ing matrix gain equation represents the diagram in Fig. 3.3. 
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(3.31) 

The matrix gain equation will be recomputed for the impaired 
aircraft. 

B. AIRCRAFT MODEL DEVELOPMENT 

A 3-axis control law model has been developed which \\rill 
operate on the stick and rudder inputs and the motion feed- 
back signals. The model outputs the control signals which 
are distributed to the aircraft actuators through the varia- 
ble gain matrix. The next step in building the F/A-18 model 
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Figure 3.3 Control Law Command Signal Distribution and Gains 
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is to develop the state variable equations for the airframe, 
actuators, and sensors. 

1 . Airframe Model 

The airframe model, obtained from the Flight Systems 
Branch at NATO, is a linearized small perturbation model for 
both the longitudinal and 1 a t e r a 1 -d i r e c t ion a 1 modes. The 
dynamic stability and control derivatives were generated from 
NATO’s F/A-18 Simulation Package. The aircraft state for the 
model is trimmed, unaccelerated 1-g flight. The airframe 
state equations are given as: 
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Airframe variable definitions and units are listed in Table 
3.3. The definitions of the stabilit}^ and control derivatives, 
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TABLE 3.3 



AIRFRAME VARIABLE DEFINITIONS 



Notation 


Variable 


Units 


u(t) 


longitudinal velocity perturbation 


f t / s 


w( t ) 


perturbed normal velocity 


f t / s 


q(t) 


perturbed pitch rate 


r a d / s 


v(t) 


perturbed lateral velocity 


f t / s 


yr(t) 


aircraft perturbed yaw rate 


r ad/ s 


rr(t) 


aircraft perturbed roll rate 


rad/s 


thed ( t ) 


aircraft perturbed pitch angle 


rad 


phi ( t ) 


aircraft perturbed roll angle 


tr 


dstx ( t ) 


collective- stabilator deflection 


tt 


d 1 e X ( t ) 


collective leading edge flap deflection 


ft 


d tex ( t ) 


collective trailing edge flap deflection 


dsty h ) 


differential stabilator deflection 


It 


d 1 e y ( t ) 


differential leading edge deflection 


tt 


dtey ( t ) 


differential trailing edge deflection 


rt 


da ( t ) 


aileron deflection 


ft 


dr(t) 


rudder deflection 


tt 


nz ( t ) 


normal acceleration 


f t / s 2 


aa(t ) 


angle of attack 


rad 


ny ( t ) 


lateral acceleration 


f t /s 2 



and associated units, which make up the matricies in Eqs. 
3.32 through 3.35 are given in Appendix E. Note that the 
units of the perturbation model are not compatible with the 
actuators or sensors. (The actuator output units are in 
degrees, and the sensor input units are in degrees/sec, de- 
grees, and G’s.) The input and output variables for the 
perturbation model were scaled in the simulation program to 
properly interface the models. Combining equations 3.32 - 
3.35 gives 
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To study reconf igurable flight controls the aircraft 
model should be capable of using the full set of control sur- 
faces available to produce the required forces and moments. 

To achieve this the control surfaces are split into indepen- 
dent right and left hand complements (i.e., right elevator, 
left elevator, etc.). The equations are then coupled so that 
a complement of control surfaces used either collectively, 
differentially, or as a single side, will produce the appro- 
priate moments. For example the stabilators deflected col- 
lectively will produce a pitching moment, deflected 
differentially will produce primarily a rolling moment, and a 
single side deflected will produce, to some degree, moments 
about all three axes. 

The F/A-18 airframe modeled in equations 3.36 and 
3.37 inherently offers control coupling through the stabila- 
tor, leading edge flap, and trailing edge flap surfaces. 
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Additional longitudinal coupling could be achieved with the 
ailerons and rudders. At the time this thesis was written 
control derivatives were not available on the longitudinal 
effects of the ailerons or rudders. (The rudder is capable 
of toe-in or flare-out and will effect the longitudinal re- 
sponse of the aircraft. This feature is normally used during 
takeoff and landing.) 

To split the control surface deflections into right 
and left hand complements, the following equations are used 
which compute the deflection inputs to Eqs. 3.36 and 3.37. 
(Also refer to Fig. 2.2 which shows the control surface posi- 
tions and corresponding positive deflections.) 

LONGITUDINAL DEFLECTIONS LATERAL-DIRECTIONAL DEFLECTIONS 
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Where r and 1 correspond to right and left surfaces. These 
equations are rewritten in the following matrix format: 
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LONG’ distribution matrix 



" d s t x" 




5 


. 5 


0 


0 


0 


0 


0 


0 


0 


O' 


d lex 


= 


0 


0 


. 5 


. 5 


0 


0 


0 


0 


0 


0 


_d t ex 




0 


0 


0 


0 


. 5 


. 5 


0 


0 


0 


0 



d s t r" 
d s 1 1 
d ler 
dlel 
d t er 
d t el 
d ar 
dal 
d r r 
drl 



(3.38) 



’LATD’ distribution matrix 



”d s t y" 




--1 


1 


0 


0 


0 


0 


0 


0 


0 


O' 


d 1 ey 




0 


0 


1 


-1 


0 


0 


0 


0 


0 


0 


d t e y 


= 


0 


0 


0 


0 


-1 


1 


0 


0 


0 


0 


da 




0 


0 


0 


0 


0 


0 


-1 


1 


0 


0 


-dr 




0 


0 


0 


0 


0 


0 


0 


0 


. 5 


. 5_ 



d s t r 
ds tl 
dler 
dlel 
"d ter 
d t el 
dar 
dal 
d r r 
drl 



(3.39) 



Note these equations are for the unimpaired aircraft only! 



Damage to one or more of the control surfaces will change the 
LONG and L ATD ma t r i c i e s directly. 



Replacing the input vectors in Eqs. 3.36 and 3.37 



with the r.h.s. of equations 3.38 and 3.39 gives the follow- 



ing modified airframe equations: 

Fm Gm 













d s t r 


u( t ) 




u(t) 






dstl 


w( t ) 


4x4 1 4x4 


w( C ) 


4x3 


3x10 


dler 


q( t ) 


Fx j 0 


q(t) 


Gx 


LONG 


dlel 


t h’e d ( t ) 


= 


C hed ( t ) 


+ 




d t er 


v(t) 


4x4 1 4x4 


v(t) 


4x5 


5'xlO 


d t el 


r ( t ) 


0 I Fyz 


yr(t) 


G y z 


LAID 


dar 


P( t) 




rr( t) 






dal 


phi( t ) 




phi ( t ) 






d r r 












drl 



(3.40) 
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Hm 








Dm 


d s t r 








u( t ) 






d s 1 1 


q ( t ) 


3x4 1 


3x4 


w( t ) 


3x3 


3x10 


d 1 e r 


n z ( t ) 


Hx 1 


0 


q(t) 


Dx 


LONG 


dlel 


aa( t ) = 






t hed ( t ) 


+ 




d ter 


y r ( t ) 


3x4 ! 


3x4 


v( t ) 


3x5 


5x10 


dtel 


rr( t) 


0 1 


Hyz 


yr(t) 


D y z 


LATD 


da r 


ny ( t ) 






rr(t) 






dal 








phi ( t ) 






d r r 














drl 



(3.41) 



The names of the individual 



matricies appear above the 



equations . 

2. Actuator Model 



Transfer functions for the flight control actuators 
were given in Reference 3 and are listed in Appendix F, 

These transfer functions are low order approximations of the 
more complicated actuator models presented in Reference 3. 

The models were used in the F/A-18 rigid body stability anal- 
ysis described in Reference 3, and approximate the frequency 
characteristics of the higher order models out to 5 hz. 

To incorporate the actuators into the modified air- 
frame model, the transfer functions are first put into state 
variable form. This procedure is outlined in Appendix F. 

The following equations represent the actuator state variable 
model : 



24x1 24x24 24x1 24x10 10x1 

Xa(t) = Fa Xa(t) + Ga Ea(t) (3.42) 

10x1 12x24 24x1 

D(t) = Ha Xa(t) 



(3.43) 



! 



where , 

Ea(t)^ = [estr(t) estl(t) eler(t) elel(t) eter(t) etel(t) 

ear(t) eal(t) err(t) erl(t)] 



and, 

D(t)^ = [dstr(t) dstl(c) dler(t) dlel(t) dter(t) dtel(t) 

dar(t) dal(t) drr(t) drl(t)] 



The input vector, Ea(t), represents the input signal in de- 
grees from the control laws, via the GAIN matrix. The output 
vector, D(t), represents the right and left control surface 
deflections in degrees. Equations 3.40 - 3.43 are combined 
to give the following equations of the airframe plus actuator 
model: 



Xp(t) Gp 



8x1 


8x8 


1 8x10 10x24 


8x1 


8x10 


Xm ( t ) 


Fm 


I Gm Ha 


Xm( t) 


0 

24x10 


24x1 


24x8 


1 24x24 


24x1 


Xa(t) 


0 


1 


Xa(t) 


Ga 



10x1 
Ea( t ) 



(3.44) 



Hp Xp( t ) 

q(t) 

nz(t) 8x1 

aa(t) 6x8 | 6x12 12x24 Xm(t) 



rr(t) 24x1 

ny(t) Xa(t) 



(3.45) 
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As before the names of the matricies appear above the equa- 
tions. The airframe plus actuator model inputs the command 
signals to the actuators, and outputs aircraft motion. Since 
the control surface deflections in the perturbation model are 
in radians, it was necessary to scale the Hm and Dm matricies 
to interface with the actuator model which outputs deflections 
in degrees. This is done in the simulation program. 

3 . Sensor Model 

In the final step of the development of the aircraft 
model, the state variable model for the aircraft sensors is 
incorporated into the airframe plus actuator model. The 
sensor transfer functions are given in Appendix G. In state 
variable form the sensor equations are given as: 

q( t ) 
nz( t ) 

11x1 11x11 11x1 11x6 aa(t) (3.46) 

X's(t) = Fs Xs(t) + Gs r(t) 

p( t) 
ny ( t ) 



eq( t ) 
enz ( t ) 

eaa(t) 6x11 11x1 

eyr(t)= Hs Xs(t) (3.47) 

e r r ( t ) 
eny ( t ) 

The sensor model inputs the aircraft motion variables in 
degrees/second, degrees, or G's, and outputs the correspond- 
ing signals to the control laws in the same units. Combin- 
ing Eqs. 3.44 - 3.47 gives the following model of airframe 
plus actuators plus sensors: 
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Fps 




X ps ( t ) 


Gps 




3 2x1 


32x32 1 


32x11 


32x1 


32x10 




X'p ( t ) 


FP 1 


0 


Xp(t) 


Gp 


10x1 
Ea( t 


11x1 


11x6 6x32 1 


11x11 


_j_ 

11x1 


11x10 


X s ( t ) 


Gs Hp 


Fp 


Xs(t) 


0 





(3.48) 





Hps 


Xps(t) 


eq ( t) 




32x1 


e n z ( t ) 


6x32 1 6x11 


Xp(t) 


eaa ( t ) = 


0 1 Hs 




eyr( t) 




11x1 


err( t) 




Xs(t) 


eny ( t) 







(3.49) 



The aircraft model inputs the actuator signals from the con- 



trol laws via the GAIN matrix and outputs the motion signals 



from the sensors which are sent to the control law .equations. 



C. ASSEMBLING THE OVERALL SYSTEM MODEL 

A mathematical model for each component in the control 
system shown in Fig. 2.3 has now been developed. Before the 
individual components of the model can be assembled, the dis- 
crete state equations for aircraft model (Eqs. 3.48 and 3.49) 
must be computed. Performing this operation the discrete 
state equations for the aircraft are given as: 



43x1 43x43 43x1 43x10 10x1 

Xps(k: + 1) = Aps Xps(k) + Bps Ea(k) (3.50) 



6x1 6x44 43x1 

^Es(k) = Hps Xps(k) 



(3.51) 



Where 



Xps(k)^ = [Xm(k) ! Xa(k) 1 Xs(k) ] 

Ea(k)'" = [esfcr(k) estl(k) eler(k) elel(k) eter(k) etel(k) 
ear(k) eal(k) err(k) erl(k)] 

Es(k)^ = [eq(K) enz(k) eaa(h) eyr(k) err(k) eny(k)] 

The Aps and Bps discrete matricies are computed as follows; 

Aps = (3.52) 

Bps = / (Is X Gps (3.53) 

Where ts represents the system sampling time. 

The GAIN matrix is now introduced to interface the dis- 
crete aircraft equations with the control law equations. Re- 
placing the input vector, Ea(k), in the aircraft equation 
with the r.h.s. of the GAIN equation (Eq. 3.31) gives the 
following : 

43x1 43x43 43x1 43x10 10x8 8x1 

Xps(k+1) = Aps Xps(k) + Bps GAIN Ec(k) (3.52) 

6x1 6x44 43x1 

Es(k) = Hps Xps(k) (3.53) 

For convenience the control law equations are repeated 
below: 

12x1 12x12 12x1 12x6 6x1 12x3 3x1 



Xc(k+1) = 


Ac 


Xc(k) + 


Bfc 


Es ( k ) + 


Be 


P(k) 


(3.29) 


8x1 

Ec(k) 


8x12 

Cc 


12x1 

Xc(k) + 


8x6 

Dfc 


6x1 

Es ( k ) + 


8x3 

Dc 


3x1 

P(k) 


(3.30) 
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Where : 



Xc(k-) ^ = [Xx(k) 1 Xz(k) ] 

EsCk)*" = [eq(h) enz(k) eaa(h) eyr(k) err(k) eny(k)] 

PCk)"^ = [px(k) py(k) pz(k)] 

Ec(k)^ = [estx(k) elex(k) etex(k) esty(k) eley(k) etey(k) 
ea(k) er(k)] 

Equations 3.52, 3.53, 3.29, and 3.30 are now combined to 
give the following matrix equation: 



43x1 


43x43 


43x10 10x8 


8x6 


6x43 1 


43x10 


10x8 


8x12 43x1 


Xps(k^l) 


A ps 


-t- Bps GAIN 


Dfc 


Hps 1 

\ 


Bps 


GAIN 


Cc Xps ( k) 


12x1 




12x6 6x43 


1 


12x12 


12x1 


Xc(k+1) 




B f c Hps 




1 




Ac 


Xc(k) 






43x10 10x8 


8x3 














Bps GAIN 


Dc 


px(k) 










+ 


12x3 




py ( k) 

pz ( k) 






(3.54) 






Be 












6x1 


6x43 


I 6x12 43x1 




6x3 







Es(k) Hps I 0 Xps(k) 0 px(k) 

= + py(k) (3.55) 

8x1 8x66x43 1 8x12 12x1 8x3pz(k) 

Ec(k) Dfc Hps I Cc Xc(k) Dc 



These equations model the dynamic response of the F/A-18 
system shown in Figure 2.3. 



D. MODELING EFFECTOR IMPAIRMENT 

Effector impairment is divided into four groups termed 
’effector impairment classes’ (EIC) [Ref. 1]. Figure 3.4, 
copied from Reference 1, defines the EIC’s and indicates 
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Figure 3.4 Effector Impairment Classes 
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their affect on aircraft stability and control. As shown in 
Fig. 3.4, effector impairment class one is unique in that 
only aircraft control is affected. The remaining impairment 
classes affect aircraft stability as well as control. A 
self repairing control system must be capable of detecting 
and classifying effector damage to compute the proper recon- 
figuration gains. 

For the case of EIC=1, one or more of the aircraft con- 
trol derivatives will be altered. In the F/A-18 model de- 
veloped above this is reflected in the LONG and LATD matricies 
(Eqs. 3.38 and 3.39). For example if the right stabilator is 
impaired the elements L0NG(1,1) and LATD(1,1) would be set to 
zero: 



LONG 



0 . 5 

0 0 

0 0 



0 0 

.5 .5 

0 0 



0 0 
0 0 

.5 .5 



0 0 
0 0 
0 0 



0 0 
0 0 
0 0 



LATD 



0 1 
0 0 
0 0 
0 0 
0 0 



0 0 

1 -1 
0 0 

0 0 

0 0 



0 0 

0 0 

-1 1 

0 0 

0 0 



0 0 

0 0 

0 0 

-1 1 

0 0 



0 0 
0 0 
0 0 
0 0 
.5 .5 



Compare the above matricies with Eqs. 3.38 and 3.39 for the 
undamaged aircraft. Note that the damaged system will now 
produce a lateral input for a given longitudinal command and 
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vice versa. For the model developed in this thesis, only 
class one effector impairments are considered. 

E. CONCLUSION 

In this chapter the mathematical models for each compo 
nent in the block diagram of Fig. 2.3 were developed. The 
individual models were then assembled to form the complete 
model of the F/A-18 dynamical system. Next the simulation 
program is developed to compute the model matricies in Eqs 
3.54 and 3.55, and compute the response of the system to 
stick and rudder inputs. 
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IV. 



PROGRAM DEVELOPMENT AND MODEL VALIDATION 



A. INTRODUCTION 

To validate the F/A-18 system model a computer program 
was written to compose the model matricies in Eqs. 3.54 and 
3.55, and compute the system response to stick and rudder 
inputs. In addition the program simulates an actuation loss 
of the right or left stabilator. The program was written in 
VS Fortran on the IBM 3033 computer at the Naval Postgraduate 
School. The program is organized to offer flexibility for 
future development and modification. 

B. PROGRAM STRUCTURE 

The program may be divided into four major operations: 

1 ) Data input 

2) Air data computations 

3) Composition of the system matricies 

4) System response computation 

The operations are performed by a series of subroutines as 
shown in the flow diagram in Fig. 4.1. In Fig. 4.1 the 
solid lines indicate control flow, and the dotted lines indi- 
cate data flow. The subroutines perform the steps outlined 
in Chapter II to compose the system matricies and compute 
the response. A functional description of each subroutine 
shown in Fig. 4.1 and associated variables is given in Appen- 
dix H. Existing subroutines at NPS were used to perform the 
required matrix manipulations [Ref. 5]. These subroutines 
are also defined in Appendix H. 
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Figure 4.1 Simulation Program Flow Diagram 
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■ Execution of the program is controlled by the 'F18' Exec 
program. The Exec defines the input and output data files as 
shown in Table 4.1. No user options are provided to define 
data input/output files. The result of each subroutine com- 
putation is written to the FA18 RESULT file, and can be 
viewed by the user after program execution. A copy of the 
Exec program is given in Appendix I. 



TABLE 4.1 

INPUT/OUTPUT FILE DEFINITIONS 



Record 

Number 


Device 


File 

Name 


Type 


Use 


01 


Disk 


FA18 


DATA 


Contains the FA18 program 
input data 


02 


Disk 


FA18 


RESULT 


Contains the result of 
the operations performed 
by each subroutine 


03 


Disk 


OPTMATD 


DATA 


Contains FA18 system 
matricies (Eqs. 3.54 and 
3.55) formated for the 
control system design 
package at NPS . 


04 


Disk 


OPTPLOT 


DATA 


Contains the system 
response data computed 
by the FA18 program 
formated for the inter- 
active plotting program 
at NPS. 



1 . Data Input 

Subroutines ^ FLITEl ^ and ^ FLITE2 ^ perform the data 
input operations. All data is read from the FA18 data file 



except the control and failure parameters which are input 
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^FLITEl’ reads in the aircraft steady state flight 
conditions and control and failure parameters. As indicated 
in Fig. 4.1 the flight conditions are read from the FA18 
DATA file. The control and failure parameters are read in 
interactively. This allows the user to conveniently run the 
program for various control inputs and control surface fail- 
ures. ^FLITE2’ reads in the basic airframe stability and 
control derivatives from the FA 18 DATA file. The stability 
and control derivatives are arranged in the matrix format 
shown in Eqs. 3.33 - 3.36, and in Appendix E. 

2 . Air Data Computations 

Prior to composing the system matricies, air data 
computations are performed by subroutine ^AIRDAT’. The sub- 
routine computes the air data inputs to the control law func- 
tions using the standard atmosphere equations [Ref. 6], The 
program does not include the logic to compute atmospheric 
conditions above the gradient (troposphere) region. There- 
fore the computations are valid only up to 36,000 feet. 

3 . Composing the S v s t em Ma t r i c i e s 

As seen in Fig. 4.1 the system matricies are com- 
posed by ten subroutine operations. The operations and 
associated equations in Chapter III are as follows: 

1) Compose the LONG and LATD matricies Eqs. 3.38 & 3.39 

2) Qompose the modified airframe matricies Eqs. 3.40 & 

3.41 

3) Compose the actuator matricies Eqs. 3.42 & 3.43 
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4) Compose the sensor matricies Eqs. 3.46 & 3.47 

5) Compose the airframe plus actuator matricies Eqs. 

3.44 & 3.45 

6) Compose the airframe plus actuator sensor matricies 
Eqs. 3.48 & 3.49 

7) Compute the discrete system matricies Eqs. 3.50 & 3.51 

8) Compose the GAIN matrix Eqs. 3.31 

9) Compose the control Law matricies Eqs. 3.29 & 3.30 

10) Compose the total system maricies Eqs. 3.54 & 3.55 

The basic component matricies are composed by first generat- 
ing a null matrix of proper dimensions, and then assigning 
the coefficient values to the proper elements. These matri- 
cies include: 

1) LONG and LATD 

2) Actuator 

3) Sensor 

4) Control law 

5) GAIN 

The coefficient values to actuator and sensor models are 
written into the program as constants; they are not contained 
on a separate data file. The coefficients in the LONG, LATD, 
control law, and GAIN matricies are first computed, then 
assigned to the appropriate matrix element. 

The LONG and LATD matricies are composed by subrou- 
tine ’LONLAT’. As explained in Section III.C these matricies 
reflect control surface damage for one of the four control 
surface impairment classes. In the present program only 
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class one effector impairments are simulated. For EIC=1 the 
appropriate control surface coefficients in the ^LONG^ and 
' LATD ^ matricies are set to zero. ^LONLAT^ contains the logic 
to impair the right or left stabilator. 'LONLAT^ also com- 
poses the unimpaired airframe control matrix, GmO. This ma- 
trix will be used to compute the impaired gain matrix in 
subroutine ^VGAIN’ . 

The control law matricies are composed by subroutine 
^LAWS\ The subroutine may be divided into four operations: 

1) Compute function gains 

2) Compute the filter coefficients 

3) Compute the matrix coefficients 

4) Assign coefficients to the control law matricies as 
in Eqs. D.l & D.2. 

The gains are computed according to the function definitions 
given in Appendix A. The functions use the air data computed 
in ^AIRDAT^, steady state AOA, and normal acceleration to 
compute the gains. All functions are programmed exactly as 
shown in Appendix A except for the following non-linear 
functions: 

Function 20 Longitudinal stick gradient 

Function 1 Lateral stick gradient 

Function 14 Rudder pedal gradient 

Function 42 RSRI non-linear gradient 

These functions were programmed using the linear terms of a 
Taylor Series Expansion about the origin. Some of the func- 
tions computed in the program are not used in the control law 
model. They are provided for user information on system per- 
formance. These functions are: 
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Function 24 
Function 25 
Function 27 
Function 29 
Function 37 
Func tion 41 
Function 112 
Function 113 



Trailing edge flap schedule 
Trailing edge flap schedule qc limit 
Leading edge flap schedule 
Leading edge flap schedule ri limit 
Nz limit on AOA feedback 
Rolling surface limit schedule 
Lateral acceleration gain 
Lateral acceleration gain 



The filter coefficients were computed using the procedures 



outlined in Appendix B to transform the analog filters into 



digital filters. The control law coefficients were computed 
using the procedures outlined in Appendix C for transforming 
the control path transfer functions into state variable form. 
The function and coefficient values are output to the FA18 
RESULT file. 



Subroutine ^VGAIN^ composes the GAIN matrix based on 
control effector impairment. For the unimpaired system the 
GAIN matrix appears as in Eq . 3.31. The subroutine is de- 

signed to implement the reconfiguration algorithm described 
in Reference 1. 



Prior to composing the total system matricies the 
discrete form of the continuous system matricies must be 
computed. This is done by subroutine ’EXPINT’ [Ref. 5]. 
’EXPINT’ computes the matrix exponential, 



F p s t s 
e 



and the integral, 



t s 

/ 

o 



ef’PS'-'sds 
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The discrete system equations are then computed as: 



Aps 



F p s t s 
e 



and 



Bps 



^ ^ F p s s 
o 



ds X Gps 



The Aps and Bps matricies are written to the OPTMATD DATA 
file. The data file is formatted for the control system 
design package at NFS. 

4 . Response Computations 

System response is computed for the recursive equa- 
tion 



55x1 55x55 55x1 55x3 3x1 

X(k+1) = AF18 X(k) + BF18 U(k) 



The response is computed for 500 data points. The AF18 and 
BF18 matricies are defined in Eq. 3.54. The state and input 
vectors are defined as: 

X(h)^ = [ Xm(k) I Xa(k) | Xs(k) | Xc(k) ] 

U(k)^ = [ px py pz ] 



Response data for all 55 states is written into the 
OPTPLOT DATA file which is formatted for the interactive 
plotting routines at NPS. The variables contained in the 
state vector which are relevant to the user for viewing are: 
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X 1 = u 
x2 = w 
X 3 = q 



X 7 = p 
x8 = phi 



x21 = dter 
x23 = d tel 
x25 - dar 
X 2 7 = dal 
X 29 = d r r 
X 3 1 = d r 1 



x5 = V 
X 6 = r 



x 4 = t h e d 



x9 = d s t r 
xl3 = dstl 
xl7 = dler 
xl9 = dlel 



The output equation (Eq, 3.55) was not programmed. 

C. PROGRAM TESTING AND MODEL VALIDATION 

To test the simulation program, four sets of runs were 
made for the following control inputs: 

1) Positive longitudinal stick, no failure 

2) Positive longitudinal stick, right stabilator failed 

3) Positive lateral stick, no failures 

4) Positive rudder pedal, no failures 

All deflections were -1 inch step inputs of 3 second duration. 
All runs were made at .6 mach/10000 feet. The sampling rate 
was set to 80 hz. (Additional runs were made at 40 hz sam- 
pling rates with no noticable difference in the output re- 
sponse.) Aircraft steady state data, including the stability 
and control derivatives, are given in Appendix J. The FA18 
RESULTS file, including the function and control matrix co- 
efficient values, is given in Appendix K. 

The model was verified for: 

1) Correct direction of motion of control surfaces and 
corresponding aircraft motion. 

2) Expected aircraft response for a right stabilator 
failure. 

3) Proper augmented aircraft motion. 
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Reference 7 contains response plots for the aircraft for simi- 
lar flight conditions. These plots were used to verify the 
model for the proper augmented motion. Some of the plots 
were reproduced for the thesis and are given in Appendix L. 

The control inputs used in Reference 7 are slightly different 
than those used for thesis model. Some comparisons however 
can still be made. All response plots for the thesis model 
are given in Appendix M. 

The aircraft response to a positive longitudinal input 
for a no failure condition is shown in Figs. M-1 - M.3. 

The motion of the right stabilator, and leading and trailing 
edge flaps compare favorably to the McDonnell model response 
shown in Fig. L.l. In the thesis model the right stabilator 
initially travels -2.5 degrees. After a .5 second transient 
period the stabilator continues to -3.6 degrees in 3 seconds. 
The corresponding McDonnell response shows an initial travel 
of -2.2 degrees. After the .5 second transient period the 
stabilator continues to -1.8 degrees in 3 seconds. The 
major difference between the two responses is the stabilator 
rate, the thesis model being slightly greater than the 
McDonnell model. (Possible explanations for the discrepancies 
in the model are given below.) 

The leading and trailing edge flaps, driven by angle of 
attacks feedback, deflect approximately 7.5 degrees in 3.4 
seconds as shown in Figs. M.l and M.2. The corresponding 
McDonnell response show the leading and trailing edge flaps 
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to deflect approximately 4 degrees in 3.5 seconds. (Note the 
initial conditions of the flaps show Fig. L.l are not zero 
as in Figs. M.l and M.2.) Aircraft pitch rate and pitch 
angle response for the thesis model are shown in Figs. M.2 
and M.3. The thesis model achieved a maximum pitch rate of 
8 degrees/ second compared to 4 degrees/second for the 
McDonnell model. Also the McDonnell model achieved a con- 
stant pitch rate in 2.5 seconds. The thesis model achieved 
a 2.5 degree pitch attitude in 5 seconds compared to 15 
degrees shown in. Fig. 1.1. The increased travel of the 
stabilator shown in Fig. M.l (compare Fig. L.l) may explain 
the discrepancies in the pitch and flap responses. 

Figure M.4 shows the response of the aircraft with the 
right stabilator failed. Note the decrease in pitch rate 
and pitch angle magnitudes. Figure M.4 shows a negative 
roll rate which is to be expected for a right stabilator 
failure. Reference 7 did not give aircraft response for 
control surface failures. 

Aircraft response to a positive lateral stick input is 
shown in Figs. M.5 - M.8. The motion of the stabilator, 
aileron, and trailing edge flap compare favorably to the 
McDonnell response shown in Fig. L.2. The aileron in the 
thesis model achieves a steady state deflection of 1.1 
degrees in .5 seconds. The McDonnell response shows a steady 
state deflection of 1.6 degrees in .5 seconds. The thesis 
model stabilator (Fig. M.5) deflects 0.45 degrees in 0.5 
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seconds, the corresponding McDonnell response shows a 0.6 
degree deflection in 0.5 seconds. The trailing edge flap 
in Fig. M.6 deflects 0.4 degrees in 0.3 seconds which com- 
pares favorably to the response shown in Fig. L.2. (Leading 
edge differential flaps are not used at the flight conditions 
tested. Function 93, which sets the gain on the leading edge 
flap path, is computed as zero.) The rudder response (Fig. 
M.6) shows an initial negative deflection which is expected 
for coordinating a right turn. The magnitude of the deflec- 
tions for the two models is approximately the same, however 
the thesis response is slightly more oscillatory than the 
McDonnell response. The aircraft roll rate, yaw rate, and 
bank angle response are shown in Figs. M.7 and M.8. Again 
only slight discrepancies exist between the two models. 

The aircraft response to a positive rudder pedal input 
is shown in Figs. M.9 and M.IO. The corresponding McDonnell 
response is shown in Fig. L.3. Comparison of the two models 
shows the shape of the responses to be approximately the 
same. However the magnitudes in the thesis model are very 
much less than the McDonnell model. For example the rudder 
achieves a maximum deflection of 0.07 degrees as shown in 
Fig. M.9. Compare this to a 1.5 degree deflection for the 
McDonnell model shown in Fig. L.3. The response of the 
differential aileron and stabilator for the thesis model 
(Figs. M.9 and M.IO) is essentially zero. This is expected 
since the rudder to rolling surface interconnect gain. 
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function 39, is computed as zero for low angle of attack. 

(See Appendix A, function 39.) It is unknown why the ailerons 
and stabilators respond as shown in Fig. L.3. Possibly the 
response shown in Fig. L.3 is based on a different flight 
control program than that used in the thesis model. 

Possible explanations for the discrepancies in the re- 
sponse between the two aircraft models are: 

1) Improper derivation of the simplified control law 
model based on the assumptions and procedures given 
in Chapters I and II. 

2) Programming discrepancies resulting in erroneous 
computations . 

3) Errors in the aircraft linear small perturbation 
mod e 1 . 

4) Differences between the modeling techniques used 

in Reference 7 and in this thesis. (The MCAIR model 
included all aerodynamic and control system non- 
linearities, as well as the effects due to digital 
time delay and quantization.) 

5) The McDonnell model response made available to the 
author may be based on a slightly different flight 
program than used in the thesis model. 

The author investigated each of the items listed above. It 
is felt the assumptions used in developing the control law 
model are correct given the available information on the 
F/A-18 control system. The simplifying assumptions used in 
developing the model may account for some of the differences. 
The program code was thoroughly reviewed and revealed no dis- 
crepancies. It is assumed however that this item remains a 
possible source of error. Two possible sources of error 
exist with respect to the aircraft perturbation model: 
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1) The coefficients in the NATO small perturbation 
model are computed with respect to the aircraft 
body axes (see Appendix E). This may effect the 
magnitude of the motion variables which are fed 
back to the control laws. 

2) The aircraft model developed in Reference 7 con- 
siders the offset position of the accelerometers 
from the aircraft center of gravity. It has not 
been determined if this effect is considered in 
the NATC perturbation model. 

Further investigation into these possibilities should be 
conducted . 



69 



IV. CONCLUSIONS AND RECOMMENDATIONS 



A mathematical model was developed which simulates the 
dynamic response of the F/A-18 aircraft. The model is 
designed to implement the reconfiguration gain matrix for 
the study of r econf igurable control systems. A program 
was written which composes the discrete time state variable 
matrices of the aircraft, and computes the response to 
stick and rudder inputs. The program also simulates the 
actuation loss of the right and left stabilators. Response 
plots of the thesis model were displayed and compared to 
the aircraft model developed in Reference 7. Possible 
sources of error were discussed, and it is recommended that 
further investigation into each of these areas be conducted. 

It is also recommended that the NATC program which com- 
putes the stability and control derivatives be acquired by 
the Aeronautical Engineering Department at the Naval Post- 
graduate School, to be made an integral part of this program. 
This would expand the ability of the program to simulate the 
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response of the aircraft for any flight condition, and various 
degrees of control surface damage. 
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APPENDIX A 



FUNCTION MATHEMATICAL DESCRIPTIONS 



FUNCTION I 

LATERAL ST I CK GRAD I ENT 



AUTO FLAP UP 
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FUNCTION 4 

ROLL RATE FEEDBACK GAIN SCHEDULE 
AUTO FLAP UP (Qc.Psl 




°c 

(p^i 



Ps 

(PSFI 




GAIN 

(VO-TS/^OLT) 
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FUNCTION 7 

LATERAL COMMAND GAIN SCHEDULE 
AUTO FLAP UP (Qc, Ps ) 





gain 

( VOLTS/VOLT 1 



62B BIAS 



75 



FUNCTION 10 



RUDDER COMMAND GAIN 
AUTO FLAR UR ( Qq ) 



GAIN 

(VOLTS/VOLT) 




Qc 

(PSF) 




1.0 BIAS 



GAIN 

I VOLTS/VOLT) 
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FUNCT ION 14 



RUDDER PEDAL GRADIENT 



Y 




Y : I YK9/x/ + YK 10) X 



MODE 


YK9 


YKIO 


YCAS AND Spin 
AND AFU 


0.00072 X (FUNCTION 10) 


0.234 X (FUNCTION 10) 


YDEL AND 


0.0017 


0.128 


SPIN 


0.00144 


0.468 


YCAS AND SPIN 
AND 7JT3 


0.00072 


0.234 
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FUNCTION 17 



RUDDER PEDAL COMMAND GAIN INCREMENT 
AUTO FLAP UP (AOA) 



GAIN 

(VOLTS/VOLT) 





DIAS 



GAIN 

(VOLTS/VOLTl 
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FUNCTION 20 

LONGITUDINAL STICK GRADIENT 

F20 : X(7.0 » 0.2 |X1 I 
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FUNCTION 22 

FADER ON THE SUPERSONIC COMPENSATION 

AUTO FLAP UP (Oc) 




(PSFl 



800 BIAS 




F22 
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FUNCTION 23 

STALL MARGIN GAIN ON PITCH 
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FUNCTION 24 



TRAILING EDGE FLAP SCHEDULE 
AUTO FLAP UP lAOA, Rj) 




ALPHA 

IDEGl 



IQC/Psl 




32.76 BIAS 



F24 

IDEGl 
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4 



FUNCTION 25 

FRAILING EDGE rLAP SCHEDULE 
AUTO FLAP UP (Qc LIMIT) 




47.636 BIAS 



(P5F 




F25 

(DEG) 
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FUNCTION 27 



LEADING EDGE FLAP SCHEDULE 
AUTO FLAP UP (ADA, R} ) 

/ 



LIMITED BY 




LIMITED BY 
FUNCTION 29L 






/ / 
/ / 




F27 

DEG 
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FUNCTION 28 



LEADING EDGE FLAP SCHEDULE 
AUTO FLAP UP (Qq LIMIT) 




oc 

(PSF) 



44.55! BIAS 




F29 

(DEG) 
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FUNCTION 29 



LEADING EDGE FLAP SCHEDULE 
auto flap up (Ri LIMIT) 




IQ^/PsI 



L I u I T 
i . 1 46 
0.7 



67.3625 ai/*S 




r 23 J 

(DEGl 



0.0 



IPST) 



I CC/P 5 ) 




n rz 9 ^ 

IDE3J 



0.0 3 : AS 



F29L = 0 DEG FOR STORES OR STORES AfC P 5 < 1200 PSF 
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FUNCTION 31 

DIFFERENTIAL TRAILING EDGE FLAP GAIN SCHEDULE 
AUTO flap up IRi , P5 ) ' 



OAIN 

(VC3LTS/VOLTI 



o.2r 




0.2247 BIAS 




GAIN 
IVOLTS/^VOLTI 



/ 
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FUNCTION 32A 



LONGITUDINAL FORWARD LOOP GAIN SCHc^DULE 

AUTO FLAP UP (Qc) 




Qc * PSF 



IPSF) 




F32A 
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1 



FUNCTION 32S 

LONGITUDINAL FORWARD LOOP GAIN SCHEDULE 
AUTO FLAP UP DthETAq£l ' 




0.13 BIAS 




K32B 



90 



1 




FUNCTION 34 



DIFFERENTIAL TRAILING EDGE FLAP 
GAIN SCHEDULE 
AUTO FLAP UP (ADA) 
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FUNCTiQN 35 

LATERAL forward LOOP GA I T4 SCHEDULE 
AUTO FLAP UP (AOA) 




AOA 

(DEG) 



0.955 BIAS 




6 DEG 



6 DEG 



LIMIT 




^0.5 


— 1 2.0 

1 


♦0.175 





gain 

(VDLTS/VOLT) 
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FUNC I I ON 35 

A 1 LERON GA I N SCHEDULE 
AUTO FLAP UP f Qq I ^5 » I 



Qr (REVERSAL POINT) SCHEDULE 



OH 

IVOLTS/VOLT) 



t .0 



0.5 



f 0 . 43 » ( 0 . ea . 

I .01 1 .01 



hinge MOMENT GAIN SCHEDULE (Rj , 

Pg ; 400 PST 
Pg * C40 PST 

Pg r 2116 PS5 




10 . 8 , 0 . 02 ) 



10,8. 0.45) 



0.5 1.0 



reversal gain SCHEDULE IQc, Gr) 
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FUNCTION 37 



NZ LIMIT ON AOA FEEDBACK 
AUTO FLAP UP (NZA) 



GAIN 

(VOLTS/VOLT) 





GAIN 

(VOLTS/VOLT) 



94 




FUNCTION 33 

Rsai gain schedule 

AUTO FLAP UP lAOA, R], Pg I 
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FUNCTION 39 



RUDDER PEDAL TO ROLL CAE INTERCONNECT 

GAIN SCHEDULE 
AUTO r L A P UP ( A 0 A 1 
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function 40 

PITCH RATH TEEDOACX GAIN SOHEXXJLE 
AUTO n_.AP UP IPs. Oc. R] I 
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FUNCTION 41 

ROLLING SURFACE LIMIT SCHEDULE 
AUTO FLAP UP (AOA.Ri) 



LJUIT 

(VOLTS) 




B.O BIAS 



I .0 BIAS 




Liwr 

(VOLT! 
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1 — 
-30 



( - IB 



INPUT 

fVOLTSl 



FUNCTION 42 

R5RI NONLINEAR GRADIENT 
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FUNCTION 45 



DIRECTIONAL FORWARD LOOP GAIN SCHEDULE 
AUTO FLAP UP (Qc, Ps > 



GAIN 

(VOLTS/VOLT) 




500 BIAS 



( P5F ) 



(P5F) 




GA I N 

I volts/volt 1 
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FUNCTION 68 



PITCH RATE FEEDBACK GAIN SCHEDULE 
AUTO FLAP UP (Qc) 




-480 BIAS 




Qc , 

IPSF 1 



LIMIT 

480 

260 



F68 



FUNCTION 90 

LATERAL ACCELERATION FEEDBACK 
GAIN SCHEDULE (R] , Ps 1 



GAIN 

VOLTS/VO-T 
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FUNCTION 93 



DIFFERENTIAL LEADING EDGE FLAP GAIN SCHEDULE 
AUTO FLAP UP (R], PS] . ^^Zaf’ 






0.69 BIAS 




gain 

I V0LTS/'*^0LT) 
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FUNCTION 100 



SPEEDBRAKE COMPENSATION INCREMENT 
AUTO FLAP UP (R ] ) 



DELTA^j 

(DEG) 




R 



4.0 BIAS 




F (00 
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FUNCTION 101 



DIFFERENTIAL STABILIZER LOAD 
auto flap up (R], 



ALLEVIATION SCHEDULE 
Ps. Nz^p) 



ga:n 

l VOLTS/VOLT) 




0.015) 

0.009) 



0.58 BIAS 




Gain 

I VOLTS/VOLT) 
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FUNCTION 96 



YAW RATE GAIN SCHEDULE 
AUTO FLAP UP (Qc, PS ) 



Qc 

(psr I 



Ps 

(PSM 



GAIN 

tVOLTS/'VOLTI 





Gain 

IVCX.TS/VO.TI 
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FUNCTION 107 



LONGITUDINAL INERTIAL GAIN SCHEDULE 
AUTO FLAP UP (Qc) 



GAIN 

I VCX.TS/VOLT ) 




°C 

ipsn 



750 BIAS B. 4X1 O' ■'bias 




GAIN 

I VOL IS /VOLT) 
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FUNCTION 108 



Qc 

tPSf) 



DIRECTIONAL INERTIAL GAIN SCHEDULE 
AUTO FLAP UP (Qc) 



GAIN 

(VOLTS/VOLT) 




0.01557 BIAS 




GAIN 

(VOLTS/VOLT) 
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FUNCTION 114 

RUDDER PEDAL COMMAND GAIN INCREMENT 
AUTO FLAP UP (R]l 






LIMIT 

A rs 








0 .8S2B6 


i--- 




♦U • jU 








♦ 0. 18 





0 . 1 ^^86 



GAIN 

I volts/volt I 
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FUNCTION 1 13 



LATERAL ACCELERATION GAIN 
AUTO FLAP UP lAOA) 




AOA. 
I DEG I 




0. 16667 




gain 

I VOLTS/VOLT) 



2.0 

BIAS 
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APPENDIX B 



A. PITCH 


DIGITAL FILTER MODELS 
RATE LEAD-LAG FILTER P2 


Using the 


•015(1+F22)S + 1 

P2(S) = (B.l) 

.015S + 1 

tustin transform ]et 

2 (Z-I) 

S = 

ts(Z+l) 



which gives 



P2(Z) = 


(P2N1) (P2N2) 

2( .015) (1+F22) + ts ts - 2( .015) (]+F22) 

[ ] z + f ] 

2( .015) + ts 2( .015) + ts 

2( .015) - ts 

Z' - r ] 

2( .015) + ts 

(B.2) 

( P2D) 



P. NORMAL ACCELERATION LAG FILTER P5 

1 

P5(S) = (B.3) 

,04S + 1 
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Using the Tustin transform: 



(P5N1) (P5N2) 

ts ts 

Z + 

.08 + ts .08 + ts 

P5(Z) = 

.08 - ts 

z - r 1 

.08 + ts 
( P5P) 



C. INTEGRATOR P9 



I 

P9(S) = 

S 



Reference 3 approximated this filter using the 
difference method. Let 



Z-1 

S = 

t sZ 



which gives 



(P9N1) (P9N2) 



P9(Z) 



tsZ + 0.0 



7 . - 1 



( P^D) 



(B.4 ) 



(P.5) 



backward 



(P. ) 
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D. LEADING EDGE AOA LAG FILTER PI 1 



1 

Pll(S) = (B.8) 

0.39S + 1 



Using the Tustin transform: 



(PllNl) (P11N2) 

ts ts 

f 1 z + 

.39(2) + ts .39(2) + ts 

Pll(Z) = (B.9) 

.39(2) - ts 

Z - 

.39(2) + ts 
( PHD) 



E. TRAILING EDGE FLAP AOA LAG FILTER PI 2 



1 

Pi2(s) = (n. 10) 

0.79S + 1 



Using the tustin transform: 



(PI 2ND (P12N2) 

t s ts 

r ] z + 

.79(2) + ts .79(2) + ts 

P12(Z) = (B. 11 ) 



.70(2) - ts 

Z 

.79(2) + ts 

( P12N3 ) 
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F. YAW RATE CANCELLER Y3 



S 

Y3(S) = 

S + 1 



Using the Tustin transform: 



Y3(Z) 



(Y3N1 ) 

2 

r ] z 

2 + ts 



(Y3N2 ) 
-2 



2 + ts 



2 - ts 

Z - 

2 + ts 

(Y3P) 



G. RSRI LEAD LAG FILTER Y5 



,75S + 1 

Y5(S) = 

. 5S + 1 



(B.l?) 



(B. 13) 



(B. 14) 



Using the Tustin transform: 



Y5(Z) 



( Y5N1 ) 

2( .75) + ts 

r 1 z + 

2{ . 5) + ts 



{Y5N2 ) 
ts - 2( .75) 
2( .5) + ts 



2( .5) - ts 

Z 

2( .5) + ts 

( Y5N3 ) 



( B . 1 5 ) 
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APPENDIX C 



SIGNAL PATH TRANSFER FUNCTIONS AND STATE SPACE MODELS 

A. PITCH RATE TO COLLECTIVE STABILATOR PATH HI ( Z ) 

From Fia. 3.1 the forward path transfer function is 

H1(Z) = [P9*P2*F68*F32a*F12 + P2* ( ( F68*F3 2a ) +F40 ) 1 (C.l) 

Let : 



A = 


= F68*F32A*F12 


B = 


= (F68*F32a)+F40 



Rewriting Eq. C.l using the filter Z-transform expressions 
for P9 and P2 gives 

A*(P9N1*Z + P9N2 ) * (P2N1*Z + P2N2 ) 

H1(Z) = 

(Z-P9D)*(Z-P2D) 

B*(P2N1*Z + P2N2) 

+ fC.2) 

Z-P2D 

Expanding Eq. C.2 



H1(Z) = 


2 

[ [A*P9N1*P2N1 + B*P2N1] Z 




[A*( P9N1*P2N2+P9N2*P2N1 ) + B* ( P2N2-P2N1 *P9D) 1 Z 




[A*P9N2*P2N2 - B*P2N2*P9Dl} / ( Z-P9D) * ( Z-P2D) 



(C.3) 
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Equation C,3 can be expressed in terms of the following 
partial fraction expansion: 



2 

bO*Z + bl*Z + b2 qstl ast2 

H1(Z) = = qst3 + + 

(Z-P9D)*(Z-P2D) (Z-P9D) (Z-P2D) 



(C.4) 



VTher e : 



qst3 = bO 

2 

astl = (bO*P9D + bl*(P9D) + b2 ) / (P9D-P2D) 

2 

qst2 = (bO*P2D + bl*(P2D) + b2) / (P2D-P9P) 

Equation C.4 can now be put in the following state variable 
form : 



xl (k+1 
x2(k+l 





estxl (k ) 




(C. 5) 



(C.6) 



The sytem of nomenclature used for the coefficients in the 
output equation (Eq. C.6) is described in Sec. III. A. 

In the remainder of the appendix only the first and last 
steps of the state variable development will be given for 
each transfer function. The equations for each of the 
coefficients are given in the simulation program listina. 
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B. NOIWAL ACCELERATION TO COLLECTIVE STABILATOR PATH H2 ( Z ) 



H2(Z) = [F12*F32a*3.5*P9*P2 + 3.5*F32a*P5] (C.l) 

x3(k+l) P9D 0 x3(k) 1 

= + nz (k ) f C . 8 ) 

x4(k+l) 0 PSD x4(k) 1 

x3 ( k ) 

estx2(k) = nzstl nzst2 + nzst3 q(k) (C.9) 

x4(k) 

C. PITCH STICK TO COLLECTIVE STABILATOR PATH H3 ( Z ) 

H3(Z) = [7 .0*F32a*F12*P9 + F20*F32a] (C.IO) 

Note: F20 defines the longitudinal stick gradient. The 

value 7.0 is a taylor series approximation 
about the origin of F20. (See F20 in App. A) 

x5(k+l) = P9D x5(k) + 1 px(k) (C.ll) 

estx3(k) = pxstl x5(k) + pxst2 px(k) (C.12) 

D. ANGLE OF ATTACK TO COLLECTIVE LEADING EDGE FLAP H4 ( Z ) 

H4(Z) = 1.328*P11 . (C.13) 

Note: The leading edge flap schedule is defined by F27. 
The value used in Eg. C.13 is the l.e. flap 
gradient. (See F27 in App. A) 

x6(k+l) = PHD x6(k) + 1 aa(k) (C.14) 

elex(k) = aalel x6(k) + aale2 aa(k) (C.15) 
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E. 



ANGLE OF ATTACK TO COLLECTIVE TRAILING EDGE FLAP H5 ( Z ) 



H5(Z) = 1.405*P12 



(C. 16) 



Note : 



The trailinq edge flap schedule 
F24. The value used in Eg, C,16 
gradient. (See F24 in App. A) 



is defined by 
is the t.e. flap 



x7(k+l) = P12D x7(k) + 1 aa(k) 

etex(k) = aatel x7(k) + aate2 aafk) 



(C.17). 
(C. 18) 



F. ROLL RATE TO DIFFERENTIAL STABILATOR PATH H6 ( Z ) 



H6(Z) = rv7*(F4+rk6t) (C.19) 
estyl(k) = rv7*(F4+rk6t) rr(k) (C.20) 
where : 



rv7 = MIN(F6*F35,F6-F101 ) 

rk6t = LIMIT(0, . 12-F4, . 12-F4) 
LL UL X 



G. LATERAL STICK TO DIFFERENTIAL STABILATOR PATH H7fZ) 



H7(Z) = 3.22*rv7*F7*(F13+F4) (C.21) 



Note: The lateral stick gradient is defined by FI. 

The value 3.22 is a taylor series approximation 
about the origin of FI. (See FI in App. A) 

esty2(k) = 3 . 22*rv7*F7* { F31+F4 ) py(k) (C.22) 
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H. RUDDER PEDAL TO DIFFERENTIAL STABILATOR PATH H8 ( Z ) 

H8(Z) = F14*l . 33*rv7*F39 (C.23) 

esty3(k) = F] 4*1 . 33*rv7*F39 pz(k) fC.24) 

I. ROLL RATE TO DIFFERENTIAL LEADING EDGE FLAP PATH H9 f Z ) 

H9(Z) = F93* (F4+rk6T) (C.25) 

eleyl(k) = F93* ( F4+rk6t ) rrfk) (C.261 

J. LATERAL STICK TO DIFFERENTIAL LEADING EDGE FLAP HIO(Z) 

HIO(Z) = 3.22*F93*F7*(F13+F4) (C.27) 

eley2(k) = 3 . 22*F93*F7* ( F31+F4 ) py(k) (C.28) 

K. ROLL RATE TO DIFFERENTIAL TRAILING EDGE FLAP PATH HI 1 f Z ) 

Hll(Z) = F31*F34*(F4+rk6t) (C.29) 

eteyl(k) = F31 *F34 * ( F4+rk6t ) rr(k) (C.30) 

L. LATERAL STICK TO DIFFERENTIAL LEADING EDGE FLAP HI 2 f Z ) 

H12(Z) = 3.22*F7*F31*F34*(F13+F4) (C.31) 

etey2(k) = 3 . 22*F7*F31 *F34* ( FI 3+F4 ) py(k) fC.32) 
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M. ROLL RATE TO AILERON PATH Hi 3 ( Z ) 

H13(Z) = F35*F36*.5*(F4+rk6t) (C.33) 

eal(k) = F35*F36*. 5*(F4+rk6t) rr(k) fC.34) 

N. LATERAL STICK TO AILERON PATH Hl4fZ) 

H14(Z) = 3 . 22*F7*F35*F36* . 5*( F13+F4 ) (C.35) 

ea2(k) = 3.22*F7*F35 py(k) (C.36) 

O. RUDDER PEDAL TO AILERON PATH H15(Z) 

H15(Z) = F14*F39*1 .33*F35*F36 fC.37) 

ea3(k) = F14*F39*1 . 33*F35*F36 pz(k) (C.38) 

P. YAW RATE TO RUDDER PATH H16(Z) 

H16(Z) = F45*F96*COS(alpha) *Y3 fC.39) 

Note: In the simulation program alpha is set to the 

steady state angle of attack in degrees. 

x8(k+l) = Y3D xP(k) + 1 yr(k) 

erl(k) = yrrl x8(k) + yrr2 yr ( k ) 
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0 . 



ROLL RATE TO RUDDER F.VrH K17(Z) 



H17(Z) = F45*F96*Y3*SIN(alpha ) 

+ F45*F38*Y5*F30*0.76*( (2*rra)+(2*rrst) ) (C.40) 



Mote: The second term on the r.h.s. of the eauation 

results from the rolling surface to rudder 
interconnect path. The vaJue, 0.76, represents 
F42, the RSRI nonlinear gradient. It is a 
taylor series approximation about the oriqin 
of F42. 



x2(k+l ) 
x3 ( k+1 ) 



Y3D 0 x2(k) 1 

+ rr (k ) 

0 Y5D x3(k) 1 



(C.41 ) 



x2 ( k ) 

er2(k) = rrrl rrr2 + rrr3 rrfk) (C.42) 

x3(k) 



R. LATERAL ACCELERATION TO RUDDER PATH H18(Z) 

(C.43 ) 
f C.44) 



H18(Z) = F45*F90 
er3(k) = nyr ny(k) 
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s. 



LATERAL STICK TO RUDDER PATH H19(Z) 



H19(Z) = F45*Y5*F38*F30*0.76*( (2*pya)+( 2*pyst) ) (C.45) 



(C.46) 
(C.47) 

T. RUDDER PEDAL TO RUDDER PATH H20(7) 



x4(k+l) = Y5D x4(k) + 1 pyfk) 



er4(k) = pyrl x4(k) + pyr2 py(k) 



H20(Z) = F45*Y5*F38*0.7&*F30*( (2*pza)+f 2*pzst) 

- F45*F14*( . 5-(F17*F114) ) (C.48) 

Note: The first term on the r.h.s. of the equation 

results from the rudder to rolling surface 
interconnect path. 

x5(k+l) = Y5D x5(k) + 1 pz(k) 

er5(k) = pzrl x5(k) + pzr2 pz(k) 



APPENDIX D 



CONTROL LAW MATRICIES 



[AC] 



xl (k + 1 ) 




P9D 1 




xl ( k ) 


x2(k+l ) 




P2D 1 




x2(k ) 


x3(k+l ) 




P9D 1 




X 3 ( k ) 


x4(k+l ) 




PSD 1 ^ 




x4(k) 


x5(k+l ) 




P9D 




X 5 ( k ) 


x6 ( k+l ) 


= 


PI ID 1 




x6 ( k ) 


x7(k+l ) 




P12d' 




x7(k) 


xa(k+l ) 




YJTT 




xB(k ) 


x9(k + l ) 




, Y3D 




x9(k ) 


xlO(k+l ) 




£_ 1 Y5D 




xl 0(k ) 


xll(k+l) 




Y5D 




X 1 1 ( k ) 


xl2(k+l ) 




1 Y5D 




xl 2(k ) 



[Bfc] 



10 0, 




10 0' 




0 1 0 1 




0 10 0 




q ( k ) 


0 0 0 1 




nz ( k ) 


0 0 1 1 




3 0 ( K ) 


0 0 1 




yr ( k ) 


1 0 0 




rr(k) 


1 0 1 0 




ny ( k ) 


0 0 10 




' 0 0 0 




10 0 0 





^[Bc] 

- 

" I 

To 0 

0 I 0 0 

~ 0 0 

I 1 0 

I 0 1 



px ( k ) 
py ( k ) 
p z ( [c ) 



(D.l) 
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APPENDIX E 



STABILITY AND CONTROL DERIVATIVE DEFINITIONS AND UNITS 

The matricies in the small perturbation model (Eqs. 3.32 
through 3.36) are defined in this appendix along with the 
dimensions of the stability and control derivatives. This 
information was obtained from the flight systems branch at 
the Naval Air Test Center. 

FX MATRIX 

1 2 3 4 

1 XU XW XO-WE -G*CTHT 

ZU ZW fZO+UB) -G*STHT 

2 

1-ZWD 1-ZWD 1-ZWD 1-ZWD 

MWD*ZU MWD*ZW MWD*fZQ+UB) MWD* { -G ) *STHT 

3 MU+ MW+ MQ+ 

1-ZWD 1-ZWD 1-ZWD 1-ZWD 

4 0 0 1 0 
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GX MATRIX 



1 



2 



3 



1 



0 



0 



0 



2 



ZDS 



1-ZWD 



ZDLF 



1-ZWD 



ZDTF 



1-ZV^D 



MWD*ZDS 

3 MDS+ 

1-ZWD 



MWD*ZDLF 

MDLF+ 

1-ZWD 



MWD*ZDTF 

MDTF+ 

1-ZWD 



4 



0 



0 



0 



HX MATRIX 



0 



0 



ZU 

1-ZWD 



ZW 

1-ZWD 



ZO 

1-ZWD 



0 



1 

UB 
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DX MATRIX 



1 

0 

ZDS 
1 -ZWD 

0 

1 

YV 

1-YVD 

NVD*YV 

NV+ 

1-YVD 

LVD*YV 

LV+ 

1-YVD 



FYZ MATRIX 



0 



ZDLF 

1-ZWD 



ZDTF 

1-ZWD 



0 



0 



2 



YR-UB 



1-YVD 



jWD* ( YR-UB ) 

NR+ 

1-YVD 



LVD*( YR-UB) 

LR+ 

1-YVD 



STHT 

rTHT 



3 



YP+V7B 



1-YVD 

NVD*(YP+WP) 

NP+ 

1-YVD 

LVD*(YP+WB) 
LP+ 

1-YVD 



4 

G*CTHT 
1 - WD 

NVD*G*CTFT 

1-YVD 

LVD*G*CTHT 

1-YVD 



0 



1 



0 



GYZ MATRIX 



1 


2 


3 


4 


5 


YPHT 


YDLF 


YPTF 


YDA 


YDR 


1-YVP 


1-YVP 


1-YVD 


] -YVP 


1 -Y^^P 


NVP*YDKT 


NVP*YPLF 


MVP* YPTF 


NVP*YDA 


NVP*YDR 



NDHT+ NDLF+ NDTF+ NPA+ NDR+ 

1-YVD 1-YVD 1-YVD 1-YVD 1-YVD 

LVD*YDHT LVD*YFLF LVD*YDTF LVD*YPA LVD*YDR 

LDHT+ LDLF+ LDTF+ LPA-^ LPP + 

1-YVP 1-Y\r 1-YVP 1-T\/P 1-YVP 

0 0 0 0 0 



HYZ MATRIX 



]. 

2 



3 



0 



0 



0 



1 



YV 

1-YVP 



YR 

1-YVP 



YP 

1-YVP 



0 
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I 



DYZ MATRIX 





1 


2 


3 


4 


5 


1 


0 


0 


0 


0 


0 


2 


0 


0 


0 


0 


0 




YDHT 


YDLF 


YDTF 


YDA 


YDR 




1-YVD 


1-YVD 


1-YVD 


1-YVD 


1-YVD 


The 


following 


abreviations were 


used in 


the above 



matricies : 

UB = BODY AXIS LONGTIIDINAL WIND 
WB = BODY AXIS VERTICAL WIND 
CTHT = COS (THETA) 

STHT = SIN( THETA) 
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D EFINITION OF DIMENSIONAL STABILITY DERIVATIVES 



HAME 


PARTIAL DERIYATIYE OFi 


WITH RESPECT TO: 


UNITS 


XU 


LWGITDOINAL FORCE 


FORWARD VELOCITY 


7^6" FT/SEC2 


xw 


LONGITUDINAL FCfiCE 


VERTICAL VELOCITY 


FT/SEC2 


XQ 


LONGITUDINAL FORCE 


PITCH RATE 


RAD/SEC2 


XV/D 


LONGITUDINAL FORCE 


VERTICAL ACCELERATION 


FT/SEC2 


XDSB 


LWGITUDINAL FORCE 


SPEED BRAKE 


RAD/SEC2 


XDTH 


LCNGITUOINAL FORCE 


THROTTLE 


PCT/SEC2 


XOS 


LONGITUDINAL FORCE 


HORIZaiTAL STABILATOR 


RAD/SEC2 


XDLF 


LONGITUDINAL FORCE 


LEADING EDGE FLAPS 


RAD/SEC2 


XDTF 


LOf^GITUDlNAL FORCE 


TRAILING EDGE FLAPS 


RAD/SEC2 


ZU 


VERTICAL FORCE 


FORWARD VELOCITY 


FT/SEC2 


ZW 


VERTICAL FORCE 


VERTICAL VELOCITY 


FT/SEC2 


ZQ 


VERTICAL FORCE 


PITCH RATE 


RAD/SEC2 


ZV^ 


VERTICAL FORCE 


VERTICAL ACCELERATION 


FT/SEC2 


ZDS8 


VERTICAL FORCE 


SPEED BRAKE 


RAD/SEC2 


ZDTH 


VERTICAL FORCE 


THROTTLE 


PCT/SEC2 


ZDS 


VERTICAL FORCE 


HORIZONTAL STABILIZER 


RAD/SEC2 


ZDLF 


VERTICAL FORCE 


LEADING EDGE FLAPS 


RAD/SEC2 


ZDTF 


VERTICAL FORCE 


TRAILING EDGE FLAPS 


RAD/SEC2 


nu 


PITCHING MOMENT 


FORWARD VELOCITY 


FT/SEC2 


MW 


PITCHING MOMENT 


VERTICAL VELOCITY 


FT/SEC2 


MQ 


PITCHING MOMENT 


PITCH RATE 


RAD/SEC2 


MWD 


PITCHING MOMENT 


VERTICAL ACCELERATION 


FT/SEC2 


MDS8 


PITCHING MOMENT 


SPEED BRAKE 


RAD/SEC2 


MDTH 


PITCHING MOMENT 


THROTTLE 


PCT/SEC2 


MDS 


PITCHING MCMENT 


HORIZaiTAL STABILIZER 


RAD/SEC2 


MDLF 


PITCHING MOMENT 


LEADING EDGE FLAPS 


RAD/SEC2 


MOTF 


PITCHING MOMENT 


TRAILING EDGE FLAPS 


RAD/SEC2 


YV 


LATERAL FORCE 


LATERAL VELOCITY 


FT/SEC2 


YVD 


LATERAL FORCE 


LATERAL ACCELERATiaJ 


FT/SEC2 


YR 


LATERAL FORCE 


YAW RATE 


RAD/SEC2 


YP 


LATERAL FCRCE 


RaL RATE 


RAD/SEC2 


YDA 


LATERAL FCRCE 


AILERW 


RAD/SEC2 


YDR 


LATERAL FORCE 


RUDDER 


RAD/SEC2 


YDLF 


LATERAL FORCE 


DIFFERENTIAL LE FLAPS 


RAD/SEC2 


YDHT 


LATERAL FORCE 


DIFFERENTIAL HGfilZ STABS 


RAD/SEC2 


YDTF 


LATERAL FORCE 


DIFFERENTIAL TE FLAPS 


RAD/SEC2 


LV 


RaLlNG MOMENT 


LATERAL VELOCITY 


FT/SEC2 


LVO 


RaLING MOMENT 


LATERAL ACCELERATION 


FT/SEC2 


LR 


RaLlNG MOMENT 


YAW RATE 


RAD/SEC2 


LP 


RaLING MOMENT 


RaL RATE 


RAD/SEC2 


LOA 


RaLlNG MOMENT 


AILERON 


RAD/SEC2 


LDR 


RaLING MaiENT 


RUDDER 


RAD/SEC2 


LDLF 


RaLING MaiENT 


DIFFERENTIAL LE FLAPS 


RAD/5EC2 


LDHT 


RaLING MOMENT 


DIFFERENTIAL HORIZ STABS 


RAD/SEC2 


LDTF 


RaLING MaiENT 


DIFFERENTIAL TE FLAPS 


RAD/SEC2 
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nFFINITION OF DIMENSIONAL STABILITY DERIVATIVES 

(Continued) 



NAME 


PARTIAL DERIVATIVE OF: 


WITH RESPECT TO: 


Mils 


NV 


YAWING MOMENT 


LATERAL VELOCITY 


FT/SEC2 


NVO 


YAWING MOMENT 


LATERAL ACCELERATION 


FT/SEC2 


NR 


YAWING MOMENT 


YAW RATE 


RAO/SEC2 


NP 


YAWING MOMENT 


RaL RATE 


RAD/SEC2 


NOA 


YAWING MOMENT 


AILERON 


RAD/SEC2 


NDR 


YAWING MOMENT 


RUDDER 


RAD/SEC2 


NDLF 


YAWING MOMENT 


DIFFERENTIAL LE FLAPS 


RAD/SEC2 


NDHT 


YAWING MOMENT 


DIFFERENTIAL HCftIZ STABS 


RAD/SEC2 


NDTF 


YAWING MOMENT 


DIFFERENTIAL TE FLAPS 


RA0/SEC2 
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APPENDIX F 



ACTUATOR TRANSFER FUNCTIONS 
A. STABILATOR 





S 2 


2(0.068) 






( ) 


+ S + 1.0 




Dst(S) 


82.9 


82.9 




Est(S) 


S 2 2(0. 


41 ) S 2 


2(0.59) 




(-( ) + 


S + 1 ] r( -) 


+ S + 1] 




36.4 36. 


4 105.3 


105.3 



The procedures outlined in Ref. 8 were used to put the 
transfer functions in state variable form. 

1) Arrange the transfer function as 

2 

b2 S + b3 S + b4 

4 3 2 

S + al S + a2 S + a3 S + a4 

b2 = 2.1377E+03 a2 = 1.6122E+04 

b3 = 2.4101E+04 a3 = 4.9559E+05 

= 1.4691E+07 a4 = 1.4691E+07 

al = 1.5410E+02 

2) Express the transfer function as a differential 
equation 

Dst(t) + al Dst(t) + a2 Dst(t) + a3 Dstft) + a4 Dst(t) 
= b2 Est(t) + b3 Estft) + b4 Est(t) 
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3) 


Comput 


,e the 


s tat 


e equati 


on coe 


f f icients 




UO 


= bO 




0 














U1 


= bl 


- 


al*U0 


= 0 












U2 


= b2 


- 


al*Ul 


- a2 


*U0 


= 2. 


1377E+03 




U3 


= b3 


- 


al*U2 


- a2 


*U1 


- a3 


*U0 = 


-3 .0532E+05 


U4 


= b4 


— 


al*U3 


- a2 


*B2 


- a3 


*B1 - 

' i 


a4*B0 = 2. 


7277E+07 


4) 


Arrange in 


i 

the followina 


I 

state 


variable 


form 










Fst 






Gst 






xl 




0 


1 


0 


0 


xl 


U1 






x2 


= 


0 


0 


1 


0 


x2 


+ U2 


Est(t) 




x3 




0 


0 


0 


1 


x3 


U3 






x4 




-a4 


-a3 


-a2 - 


al 


x4 


U4 







Hst 

xl 

Dst(t) =1000 x2 

x3 

x4 



Note input and output for all actuator models is in 
degrees. The remaining state variable models are developed 
as above. 
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B. LEADING EDGE FLAP 



Dle(S) 

Ele(S) 



1 .0 



S S 

( + 1 . 0 ) ( + 1 , 0 ) 

26.9 82.9 



xl(t) xl(t) 

= Fie + Gle Ele(t) 

X 2 ( t ) X 2 ( t ) 



Dle(t) 



Hie 



xl(t) 
x2 ( t ) 



C. TRAILING EDGE FLAP 



Dte(S) 


1.0 


Ete(S) 


S 2 2f0.71) 




( ) + 




35.0 35.0 



xl(t) 


xl(t) 




= Fte + Gte Ete(t) 


x2 ( t ) 


x2 ( t ) 



xl(t) 

Dte(t) = Hte 

x2(t) 
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i 



D. AILERON 



Da(S) 

Ea(S) 

xl(t) 
x2 ( t ) 

Da(t) 



E. RUDDER 

Dr(S) 
Er (S) 



xl (t) 
x2 ( t) 

Dr(t) 



1.0 

S 2 

( ) 

75.0 



xl(t) 

= Fa + 

x2 ( t ^ 



xl(t) 

= Ha 

x2 ( t ) 



1 

S 2 

( ) 

72.1 



xl ( t ) 

= Fr 

x2(t) 



xl ( t ) 

= Hr 

x2(t) 



2(0.59) 

S + 1 .0 

75 . 0 



Ga Ea ( t ) 



2(0.69) 

+ S + ].0 

72.1 



+ Gr Er(t) 
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The actuator state variable model matricies in Eas. 3.43 
& 3.44 are arranged as follows; 

FA MATRIX 



Fst 

Fst 

Fie 0 

Fie 

Fte 

Fte 

0 Fa 

Fa 

Fr 

Fr 



GA MATRIX 



Gst 

Gs t 

Gle 0 

Gle 

Gte 

Gte 

Gdi 

0 Ga 

Gr 

Gr 



HA MATRIX 



Hst 

Hst 

Hie 

Hte 0 

Hte 

Hie 

Hie 

0 Ha 

Ha 

Hr 

Hr 
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APPENDIX G 



AIRCRAFT SENSOR TRANSFER FUNCTIONS 



A. RATE GYRO TRANSFER FUNCTION 



S 

+ 1.0 

Y(S) 131.7 

E(S) S S 

r + 1.03 C + 1.01 

167.8 461.7 • 



B. ACCELEROMETER TRANSFER FUNCTION 



S 

+ 1.0 

Y(s) 235.8 

E(S) S 2 2(.96) 

( ) + S + 1 .0 

395.3 395.3 



C. ANGLE OF ATTACK SENSOR TRANSFER FUNCTION 



Y(S) 1.0 



E(S) S 

+ 1.0 

14.0 



Note: Transfer functions input degrees/second (or degrees) 
and output degrees/second (or degrees). The state variable 
model for the transfer functions in Eqs . 3.46 & 3.47 are 
developed using the procedures outlined in appendix F. 
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APPENDIX H 





SIMULATION PROGRA.M SUBROUTINES 


Subroutine : 


FLITEl 


Description ; 


: Reads in steady state flight conditions, stick 

and rudder commands, and control surface 
failure parameters. 



Calling sequence: CALL FLITEl (mach, alt , alpha, nz , ncont , nst , 

nstp, amp, rstf, Istf) 

Input arguements: None 

Output arguements ; 



mach 


steady state mach number 


alt 


Steady state altitude in feet 


alpha 


Steady state AOA in degrees 


nz 


Steady state normal acceleration. Normally 1.0 


t s 


Sam.pl ing time in seconds 


ncont 


Control number: l=Longi tud inal stick 

2=Lateral stick 
3=Rudder deflection 


nst, nstp 


Control start time, and control stop time. 

Nst and nstp are actually the iteration number: 
Start time ( secs . )=nst*ts 
Stop time f secs . ) =nstp*ts 


amp 


Control amplitude given as deflection in inches 


rstf, Istf 


Right and left stabilator failure parameters 
respectfully: No failure = 0 

Failure = 1 
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Subroutine 



FLITE2 



Operation: Reads in basic airframe matricies from the 

'F/A-18' data file. 

Calling arguments: CALL FLITE2 ( Fx , Gx , Hx , Dx , Fyz , Gyz , Hyz , Dyz , 

Nfx,Ngx,Nhx,Ndx,Nfyz,Ngyz, Nhyz , Ndyz ) 



Input arguments: None 



Output arguments: 

Fx, Gx, Hx, Dx Longitudinal state matricies 



Fyz, Gyz, Hyz, Dyz Lateral-directional state matricies 



Nfx, Ngx, Nhx , Ndx Two-dimensional vectors giving the 

Nfyz , Ngyz , Nhyz , Ndyz number of rows and columns of the 

respective m.atricies. Example, 
Nfx(l) = Number rows in Fx 
Nfx(2) = Number columns in Fx 
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Subroutine : 



LONLAT 



Description : 


Generates the LONG and LATD matricies 
( Eqs . 3.38 and 3.39) based on the control 
surface failure parameters. 

Computes the unimpaired input matrix, GmO, to 
the modified ai frame eauations. 'Hie GmO matrix 
v;ill be used in the reconfiguration algorithm. 



Calling Arguments: CALL LONLAT( rstf , 1st f , i fai 1 , Gx , Gyz , LONG , 

LATD , GmO , Nlong , Nlatfi , NgmO , Max , Mgyz ) 

Input arguments: 



rstf, Istf 


Right and left stabilator failure parameters 
respectfully. (See subroutine FLITEl ) 


i f ai 1 


Failure flag: No failure = 0 

Failure = 1 


Gx, Gyz 


Longitudinal and lateral-directional input 
matricies . 


Ngx, Ngyz 


Two dimensional vectors giving the number of 
rov/s and columns of the G.x, and, Gyz 
matricies . 



Output arguments: 



LONG , LATD 


Matricies described in Section III.B. which 
split the control surface deflections into 
right and left hand sides. The LONG and LATD 
matricies reflect control surface failure or 
damage . 


GmO 


The unimpaired input matrix to the modified 
airframe equations ( Eq . 3.40). GmO is 
composed using the unimpaired LONG and LATD 
matricies. It is subsequently used in the 
reconfiguration algorithm. 


Nlong, Nlatd 
NgmO 


Two dimensional vectors giving the number of 
rows and columns in the respective matricies. 
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Subroutine: AIRDAT 



Description: Computes the air data inputs to the control 

laws based on standard atmosphere conditions. 

Calling arguments: CALL AIRDATf mach, alt , temp, rho, ps i , a , 

ri) 

Input arguements: 

mach Steady state mach number 

alt Steady state altitude in feet 

Output arguements: 

temp Standard atmosphere temperature in degrees 

rankine 

rho Standard atmosphere density in 

lb-sec**2 / ft**4 

psi Standard atmosphere pressure in lb/ft**2 

a Sonic velocity in ft/sec 

qc Dynamic pressure in lb/ft**2 

ri Pressure rat io=gc /psi 
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Subroutine: MODEQ 



Description: Composes the modified airfram.e matricies (Eos. 

3.40 and 3.41) 



Calling Arguments: CALL MODEQ( Ex , Gx , Hx , Dx., Fyz , Gyz , Hyz , Dyz , 

LONG, LATD, Fm , Gm, Hm , Dm , Nf x , Ngx , Nhx , Nd x , 
Nfyz , Ngyz , Nhyz , Ndyz ) 



Input arguments: 

Longidudinal and lateral-directional 
basic airframe matrici.es (Eos. 3.32 
- 3.35) Read in from subroutine 
FLIGHTl 

LONG, LATD LONG and LATD matricies (Eos. 3.38 

and 3.39) Generated in subroutine 
LONLAT 



Fx, Gx , Hx , Dx 
Fyz, Gyz, Hyz, Dyz 



Nfx, Ngx, Nhx, Ndx 
Nfyz, Ngyz, Nhyz, Ndyz 
Nlong, Nlatd, 



Two dimensional vectors giving the 
number of rows and columns in the 
respective matrix. 



Output arguements: 



Fm, 


Gm, Hm, 


Dm 


Modified basic 
( Eqs . 3.40 and 


airframe matricies 
3.41 ) 


Nfm, 


Ngm , 


Nhm, Ndm 


Row and column vectors as described 
as described above 
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Subroutine : 



SENSOR 



Description; Composes the sensor matricies (Eqs. 3.47 and 
3.46) 

Calling arguments: CALL SENSOR( Fs , Gs , Hs , Nf s , Ngs , Nhs ) 

Input arguments; None 
Output arguements: 

Fs, Gs, Hs Sensor matricies 

Nfs, Ngs, Nhs Two dimensional row and column vectors 



Subroutine: ACTU 

Description: Composes the actuator matricies (Eqs. 3.42 and 

3.43) 

Calling arguements: CALL ACTU( Fa , Ga , Ha , Nfa , Nga , Nlia ) 

Input arguments: None 

Output arguments; 

Fa, Ga, Ha Actuator matricies 

Nfa, Nga, Nha Two dimensional Row and column vectors 
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Subroutine ; 



PLANT 



Description: Composes the basic airframe plus actuator 

matricies (Eqs. 3.44 and 3.45) 



Calling arguments : 


CALL PLANT ( Fm , Gm , Hm , Dm , Fa , Ga , Ha , Fp , Gp , Ho , 
Nf m, Ngm, Nhm, Ndm , Nfa , Nga , Nha , Nfp, Ngp, Nhp ) 



Input arguments: 



Fm, Gm, Hm, Dm 


Modified basic airframe matricies 
(Eqs. 3.44 and 3.45) Generated in 
subroutine MODEQ. 


Fa , Ga , Ha 


Actuator matricies (Eqs. 3.42 and 3.43) 
Generated in subroutine ACTU. 


Nfm, Ngm, Nhm, Ndm 
Nfa, Nga, Nha 


Two dimensional row and column vectors 



Output arguements : 



Fp , Gp , Hp 


Airframe plus actuator matricies 


Nfp, Ngp, Nhp 


Row and column vectors 
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Subroutine: LAWS 



Description; Computes the function gains (App. A) 

Computes the coefficients in the control law 
matr icies 

Composes the control law matricies 

Calling arguments; CALL LAWS ( alpha , nz , psi , qc , r i , ts , Ac , Bfc , 

Be , Cc , Df c , Dc , Nac , Nbf c , Nbc , Ncc , Nd f c , Ndc ) 

Input arguments; 

alpha Steady state angle of attack in degrees 

nz Steady state normal acceleration =1.0 

psi standard atmosphere pressure at steady state 

altitude in lb/ft**2 

qc dynamic pressure in lb/ft**2 

ri pressure ratio 

ts sampling time in seconds 

Output arguments; 

Ac, Bfc, Be Control lav/ matricies (Eqs. 3.29 and 3.30) 

Cc, Dfc, Dc 

Nac, Nbfc, Nbc Row and column vectors 
Ncc, Ndfc, Ndc 
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Subroutine : 



VGA IN 



Description: Computes the variable gain matrix based on 

control surface im.pairment. The subroutine 
is disigned to use reconfiguration algorithm 
introduced in Ref. 1. 



Calling arguments: CALL VGAIN( i fa i 1 , i f ix , Gm, GmO , GAIN, Ngm. , 

NgmO, Ngain) 



Input arguments: 



i fai 1 



i f ix 



Gm 



Control surface failure flag set in subroutine 
■LONLAT; No failure = 0“ 

Failure = 1 

Fix parameter. The GAIN matrix will be computed 
fo the impaired aircraft if ifix = 1 and 
ifail = 1. Ifix is read interactively by sub- 
routine FLITEl . 

Airframe control matrix. From subroutine MODEO. 



GmO 



Ngm, NgmO 



Unimpaired airframe control matrix. Generated 
in subroutine LOMLAT. Note if none of the 
control surfaces are damaged or failed, Gm=GmO . 

Row and column vectors 



Output Arguements: 

GAIN Variable gain matrix described in section 

III. A. 

Ngain Row and column vector 
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Subroutine: EXPINT [Ref. 5] 



Description: Computes both the matrix exponential- 

Fps*ts 

Aps = e 

and the integral 

Fps*s 

^ e ds 

O 

Reference 5 gives the method used to perform 
the computations. 

Calling arguements: CALL EXPINTf Fps , Nfps , Aps , Naps , Duml , 

Nduml , ts , iop, Dum2 ) 



Input arguments: 

Fps Airframe plus actuator plus sensor transfer 

matrix. Generated in ARCRFT 

Nfps Two dimensional row and column vector 

ts Samplina time 

iop Print parameter: 0 = Do not print results 

Fns*ts 

Otherwise print Fps, Ts , e 
Dum2 Vector of working space for computations 



Output arguments: 

Aps Matrix exponential described above 

Naps Two dimensional row and column vector 

Duml 'The integral of the matrix exponential. To be used 

in computing the discrete control matrix Bps. 

43x10 43x43 43x10 

Bps = Duml X Gps 
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The following subroutines perform the matrix operations 
necessary to com.pose the A/F-18 system matricies. 

Subroutine: MULT [Ref. 5] 

Description: Performs the matrix multiplications. 



Calling arguments: CALL .MULT( A, NA, B, NB, C, NC) 

Input arguments: 

A, B Matricies to be m>ultiplied 

Na, Nb Two dimensional vector giving the number of rows 

and columns in A, and B. 

Output arguments : 

C Product of A and B, C = AB 

Nc Row and column vector 



Subroutine ADD and SUBT are called just as MULT. 



Subroutine: NULL [Ref. 5] 



Description: Generates a null matrix 

Calling arguments: CALL NULL(A,NA) 



Input arguments: 

Na Two dimensional vevtor giving the number of rows 

and columns of the null matrix, 

Na ( 1 ) = Number of rows 
Na ( 2 ) = Number of columns 



Output arguments: 

A Null matrix having dimension Na ( 1 ) x Na ( 2 ) 



Subroutine: JUXTR [Ref. 5j 

Description: Juxtaposes two matricies by row 

Calling arguements: CALL JUXTR( A, Na , B, Nb, C, Nc ) 



Input arguements: 

A, B The two matricies to be juxtaposed by row 

Na , Nb Row and column vectors 



Output arguments: 



C 

Nc 



The juxtaposition of A and B C 
Row and column vector 



A 

B 
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Subroutine: JUXTC [Ref. 5] 



Description: Juxtaposes two matricies by column 

Calling arguments: CALL JUXTR( A, Na , B, Nb, C, Nc ) 

Input arguements: 

A, B The two matricies to be juxtaposed by column 

Na, Nb Row and column vectors 

Output arguments: 

C The juxtaposition of A and B C = A I B 

Nc Row and column vector 



Subroutines: ADD and NULT 

The calling sequences for these subroutines are exactly as 
described for those matrix operations above. 
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Subroutine: OUTPUT 



Description: Outputs a matrix with rows and columns 

numbered . 



Calling Arguments: CALL OUTPUT( A, Na { 1 ) , Ma ( 2 ) , ' A ') 



Input arguments: 

A Matrix to be output 

Na(l), Na(2) Number of rows and number of columns 

in A respectfully 

'A ' Matrix name. Must be placed in quotes 

and occupy four spaces. 

Output arguments: 

The A matrix is written to the FA18 RESULTS file 
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APPENDIX I. 



F/A-18 EXEC PROGPJ^ 



&TRACE OFF 
CLRSCRN 

GLOBAL TXTLIB VALTLIB VFORTLIB IMSLDP NONIMSL 

GLOBAL LOADLIB VFLODLIB 

FILEDEF 06 TERMINAL 

FILEDEF 01 DISK F/A-18 DATA A1 

FILEDEF 02 DISK F/A-18 RESULTS ( LRECL 132 

FILEDEF 03 DISK OPTMATD DATA A1 

FILEDEF 04 DISK OPTPLOT DATA A1 

&TYPE LOADING FA18 SIMULATION 

LOAD FA18 ORACLSl (START 

&EXIT 



CMSLIB 
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APPENDIX J 



FLIGHT CONDITIONS AND STABILITY AND CONTROL DERIVATIVES 



A. FLIGHT CONDITIONS 



MACH = 0.6 


ALT = 


10000.0 


VTRUE = 646.42 


ALPHA = 2.6184 


THRUST 


= 3531.5 


GAYA =0.0 


PCTHLC = 20.496 


PCTHRC 


= 20.496 


XIXX = 21227.0 


DLHTD = 0.68914 


DRHDT 


= 0.68914 


XIXZ = -1827.0 


DLAD =0.0 


DRAD = 


o 

o 




DLEFL = 3.0968 


DLEFR 


= 3.0968 




DTEFL = 3.2647 


DTEFR 


= 3.2647 




DRUDL =0.0 


DRUDR 


= 0.0 




FLARUD =0.0 


DELSB 


= 0.0 




CG = 0.2206 


WAIT = 


32550 




XIYY= 0.1220E+06 


XIZZ = 


0. 13958E+06 





B. LONGITUDINAL DERIVATIVES 



XU = 


-0. 13257E-01 


ZU = 


-0.73337E-01 


MU = -0.12988E-04 


XW = 


0. 71265E-01 


ZW = 


-1 .1526 


MW = -0.11331E-01 


XQ = 


0.32650 


ZO = 


-5.6525 


MO = -0.59346 


X^'JD = 


= 0.39729E-03 


ZWD = 


: -0.6917 


MWD = -0.34049E-03 


XDSB 


= 0.0 


ZDSB 


= 0.0 


MDSB =0.0 


XDTH 


= 0.14257 


ZDTH 


= 0.0 


MDTH =0.0 


XDS = 


^ -0.43034E-03 


ZDS = 


^ -0.19801E-01 


MDS = -0.24151E-02 


XDLF 


= -0.412942E-03 


ZDLF 


= 0.64114E-02 


MDLF = -0.5230E-03 


XDTF 


= 0.17674E-03 


ZDTF 


= -0.45250E-01 


MDTF = 0.32773E-01 



C. LATERAL- DIRECTIONAL DERIVATIVES 



YV = -0.2437 


LV = 


-0.29125E-01 


NV = 0.69481E-02 


YVD =0.0 


LVD = 


0.0 


NVD =0.0 


YR = 0.77567 


LR = 


0.70076 


NR = -0.21705 


YP = 0.67478E-02 


LP = 


-3.1241 


NP = -0.15254E-01 


YDA = 0.86683E-03 


LDA = 


-0.43808E-02 


NDA = 0.77305E-04 


YDR = 0.52918E-02 


LDR = 


0.72058E-03 


NDR = -0.46503E-03 


YDLF =0.0 


LDLF 


= 0.0 


NDLF =0.0 


YDHT = -0.21027E-02 


LDHT 


= 0.39478E-02 


NDHT = 0.57361E-04 


YDTF =0.0 


LDTF 


= 0.24948 


NDTF = -0.32655E-02 
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Eimm. CQMDITIQN 



VARIABLE 


DEFINITION 


LMIS 


XMACH 


MACH NUMBER 


ND 


ALT 


ALTITUDE. M3L 


FT 


VTRUE 


TRUE AIRSPEED 


FT/SEC 


ALFA 


ANGLE Of ATTACK 


DEG 


THRUST 


TOTAL ENGINE THRUST 


LB 


GAMA 


FLIGHT PATH ANGLE 


RAD 


PCTHLC 


LEFT THROTTLE POSITION CCMMAND (lOOS'MAX AB) 


PERCENT 


PCTHRC 


RIGHT THROTTLE POSITICtI COMMAND! 1 00SH1AX AB) 


PERCENT 


DLHTD 


LEFT HORIZOI^TAL STAB DEFLECTION (<■ TE DOV.7J) 


DEG 


DRHTD 


RIGHT HCRIZOfJTAL STAB, DEFLECTION (+TE DOWN) 


DEG 


DLAD 


LEFT AlLEROfJ DEFLECTlOf7 ( + TE DOWN) 


DEG 


DRAD 


RIGHT AILERON DEFLECTlOf) (+TE DOWN) 


DEG 


DLEFL 


LEFT LEADING EDGE FLAP DEFL (HE DOWN) 


DEG 


DLEFR 


RIGHT LEADING EDGE FLAP DEFL. (HE DOWN) 


DEG 


DTEFL 


LEFT TRAILING EDGE FLAP DEFL. (+TE DOWN) 


DEG 


DTEFR 


RIGHT TRAILING EDGE FLAP DEFL (+TE DOWN) 


DEG 


DRUDL 


LEFT RUDDER DEFLECTlOfJ ( + TE LEFT) 


DEG 


DRUDR 


RIGHT RUDDER DEFLECTION (+TE LEFT) 


DEG 


FLARUD 


FLARED RUDDER ( PA/TO ONLY ) 


DEG 


DELSB 


SPEEDBRAKE DEFLECTION 


DEG 


C6 


LONG. CENTER OF GRAVITY LOCATIOt-l IN PCT MAC 


PERCENT 


WAIT 


AIRCRAFT WEIGHT 


LB 


XIXX 


MOMENT Cf INERTIA ABOUT X-BOOY AXIS 


SLU6-FT2 


XIYY 


Ma-IENT Of INERTIA ABCOT Y-BOOY AXIS 


SLUG-fT2 


XIZ2 


MOMENT Of INERTIA ABOUT Z-600Y AXIS 


SLUG-FT2 


XIXZ 


XZ PLANE PRODUCT Cf INERTIA 


SLUG-FT2 
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AIRCRAFT FLIGHT CONDITIONS AND SAMPLING TIME 



APPENDIX K 



F/A-18 RESULTS FILE 
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MCDONNELL DOUGLAS MODEL RESPONSE PLOTS 
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Figure L.l Longitudinal Stick Response Mach 0.6/10,000 Feet 
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Figure L.2 Lateral Stick Response Mach 0.6/10,000 Feet 
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Figure L.2 Continued 
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Figure L.3 Rudder Response Mach 0.6/10,000 Feet 
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Figure L.3 Continued 
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Figure 



-10 RESPONSE TO 3 SEC 1 IN LONG PULSE 
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F-18 RESPONSE TO 3 SEC 1 IN LONG PULSE 
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F-18 RESPONSE TO 1.0 IN 3 SEC. 
LONGITUDINAL PULSE WITH 
RIGHT STABILATOR CENTERED AND IX 
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F-18 RESPONSE TO 1 SEC 0.5 IN LAT PULSE 
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Figure M 




-18 RESPONSE TO 1 SEC 0^ IN LAT PULSE 
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F-18 RESPONSE TO 1 SEC 0.5 IN LAT PULSE 
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F-10 RESPONSE TO 1 SEC 0.5 IN LAT PULSE 
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F-18 RESPONSE TO 1 SEC 0^ IN DIR PULSE 
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Figure M . 9 




F-10 RESPONSE TO 1 SEC 0.5 IN DIR PULSE 
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